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STUDY OF LOAD ALLEVIATION AND MODE SUPPRESSION 
(LAMS) ON THE YF-12A AIRPLANE 

By Lester D. Edinger, Frederick L. Schenk, and 
Alan R. Curtis 

Lockheed ADP, Burbank, California, and 
Honeywell GAPD, Minneapolis, Minnesota 


SUMMARY 


The intent of the study was to evaluate the potentials and capability for 
implementing a LAMS (load alleviation and mode suppression) system on the 
YF-12A for the purpose of flight research. This research would be a con- 
tinuation of work performed in LAMS technology on previous programs such 
as the NASA/Air Force XB-70 and the Air Force B-52 flight test programs 
(refs. 1 and 2). The nature of the research is to minimize the design risk in 
application of LAMS to future aircraft. 

The results of the study show that the YF-12A would be a suitable test 
bed for continuing development of LAMS technology. This was demonstrated 
by defining five candidate LAMS systems and analytically evaluating them with 
regard to performance and mechanization. Each of the five systems used a 
different combination of force producers. A small canard vane or a mass- 
reaction device mounted near the cockpit were considered as possible LAMS 
force producers, together with the existing inboard and outboard elevons. It 
was concluded that a combination of canard vane and outboard elevons would 
provide the most effective system for the YF-12A. 


INTRODUCTION 


LAMS (load alleviation and mode suppression) is a flight control techno- 
logy directed toward improving the efficiency of flight. It uses the principles 
of feedback control to regulate the loads on the flexible aircraft induced by 
either flight through turbulent air or by abrupt pilot command inputs. It 
extends the techniques used in conventional flight controls to include stability 
augmentation of the aircraft's significant aeroelastic modes. 

Most of what is considered today as LAMS technology has evolved from 
a NASA/Air Force XB-70 program and an Air Force B-52 program. These 
programs are described in refs. 3 and 4, respectively. Prior to these pro- 
grams, several smaller studies on airplanes were performed, but most of the 
prior work in load alleviation and structural mode control was confined to 
large flexible boosters such as the Saturn 5. Although the potential of this 
technology had been established for boosters, the benefits possible for the 
airplane had not been determined. However, experience with structural 



flexibility on existing aircraft such as the B-47, B-52, and XB-70, together 
with extrapolation to potentially more serious problems on vehicles such as 
the SST, C-5A, 747, and the B-1A, pointed out a need for development of 
LAMS technology. 

The objective of the XB-70 program was to evaluate the use of feedback 
control techniques to suppress the aircraft's aeroelastic contribution to local 
accelerations. This comprehensive study led to the design and flight test 
evaluation on the XB-70 of a system referred to as ILAF (identical location 
of accelerometer and force). Results of the flight test are described in ref. 2. 
Although the XB-70 was retired before the tests could be completed, limited 
results obtained at two flight conditions indicated the feasibility of active 
structural mode control. During the flight tests, a small canard vane 
(approximately four sq ft in area) was used as a means for exciting the struc- 
tural modes to test the system. The demonstrated effectiveness of the vane 
for exciting the modes led to the conclusion it could also be used as an excel- 
lent force producer for controlling the modes. However, termination of the 
XB-70 program precluded investigation of such a system. 

The B-52 program, which ran almost concurrently with the XB-70 pro- 
gram, was established to demonstrate through flight test the use of LAMS 
techniques to reduce aircraft fatigue damage rates in air turbulence. This 
represented a formidable problem, as it was necessary to evaluate stress 
at several vehicle stations and to use several force producers to achieve the 
desired control. This compounded the synthesis task; hence, the use of 
modern control theory was introduced as a design tool. The B-52 LAMS 
system was designed primarily at one flight condition but its performance 
was evaluated at three discrete conditions. The test results (ref. 1) sub- 
stantiated the performance predicted by the analytical studies. Although the 
XB-70 and the B-52 programs served to establish the potential of LAMS 
technology, additional development is required before LAMS can be applied 
without significant design risk. Examples of areas requiring further develop- 
ment include evaluation of the effects of LAMS on flutter margins, effective 
use of auxiliary force producers, and improvements in LAMS sensing techni- 
ques. Continued development is also expected to lead to new applications for 
LAMS, such as augmentation of flutter margins and maneuver load control. 

The intent of the LAMS study on the YF-12 is to continue the development 
of LAMS technology. With the retirement of the XB-70 as a test aircraft, the 
YF-12A is the most promising candidate available to NASA for a LAMS design 
base. TheYF-12A, like the XB-70, is a low -load -factor, low-aspect-ratio, 
supersonic aircraft; hence, it exhibits aeroelastic properties characteristic 
of contemporary transport aircraft. The YF-12A itself does not require 
LAMS technology, and it is not well suited for evaluating rigid-body load alle- 
viation techniques such as direct lift control. However, because of its aero- 
elastic properties, it is well suited for evaluating structural mode suppression 
techniques. With the present state of the art of load alleviation and mode 
suppression, use of the latter techniques constitutes most of the design risk. 
Hence, evaluation of mode suppression techniques on the YF-12A could result 
in eliminating most of the design risk associated with LAMS technology and, 
thus, lead to acceptance by the airframe manufacturers. 
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SYMBOLS 


A/C 

BL 

dB 

F/C 

ft/ sec 

FS 

WS 

h 

i 

K G(s) 
c 

LAMS 

M 

n 

z 

PSD 

q 

q c 

rms 

s 

SAS 

w 

W g 

r 

z 


aircraft 
butt line 
decibel 

flight condition 
feet per second ' 
fuselage station 
wing station 
altitude 

yrr 

transfer function 

load alleviation and mode suppression 
Mach number 
load factor 

power spectral density 

2 

dynamic pressure, lb/ ft 

p 

pitot differential pressure, lb / ft 
root mean square 
Laplace operator 
stability augmentation system 
vertical velocity, in. /sec 
wind state 
desired attenuation 
vertical displacement, in. 
inboard elevon deflection, rad 



SYMBOLS (concluded) 



x 


£ 

0 

0 

a 

u 


outboard elevon deflection, rad 

canard deflection, rad 

reaction mass-actuator command, in. 

mass displacement, in. 

damping ratio 

pitch angle, rad 

rigid-body pitch rate, rad/ sec 

standard deviation 

frequency, rad/ sec 

STUDY APPROACH 


The objective of the study was to establish the effectiveness of using 
feedback control techniques to suppress the aircraft's structural modes. To 
meet this objective, it was necessary to synthesize LAMS controllers and 
analytically demonstrate their potential mode suppression performance. The 
mode suppression performance was evaluated in terms of the percentage 
reduction in root-mean-square (rms) normal acceleration at selected stations 
on the aircraft. The acceleration response was evaluated for both a gust 
input disturbance and input disturbances introduced through excitation of a 
control surface. 

The analysis was restricted to the longitudinal axis of the airplane. 
Vehicle data were provided for four flight conditions. Detailed system syn- 
theses were performed at one condition, but performance and stability were 
evaluated at all four flight conditions. 

Several combinations of force producers were evaluated in the study to 
determine the most effective set for mode suppression. These combinations 
included not only the existing force producers but an auxiliary force producer 
as well. The existing force producers consist of inboard and outboard ele- 
vons. The candidates considered for the auxiliary force producer were a 
small canard vane or a mass-reaction device, either one being located near 
the nose of the aircraft. Since neither one of these auxiliary force producers 
presently exist on the YF-12A aircraft, feasibility of implementation had to be 
considered. 

Since up to three force producers could be used in a system and since 
performance was to be evaluated at several aircraft stations, the synthesis 
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was a multiple -input/ multiple-output control problem. Quadratic optimiza- 
tion techniques from modern control theory are especially suited for control 
problems of this type. Hence, both conventional frequency response analysis 
techniques and quadratic optimization techniques were used in the study to 
define the control systems. 

Once analytical block diagrams of candidate controllers were defined, 
studies were made to determine system mechanizations. Mechanization 
studies included primarily determination of potential feedback sensor location 
problems and required aircraft modifications 


AIRCRAFT DESCRIPTION AND ANALYTICAL REPRESENTATION 

Aircraft Description 

The YF-12 aircraft is an advanced two-engine, supersonic interceptor 
with a maximum gross weight of approximately 124 000 lb. The aircraft con- 
figuration is shown in Figure 1. As Figure 1 indicates, the aircraft has a 
delta wing with a chine construction on the forebody. The chine gives the 
vehicle aerodynamic properties similar to those of a double delta configura- 
tion. The two engines are located at about mid-wing on each side, dividing 
the wing into inboard and outboard sections. Each of the four sections has an 
elevon surface which is used for pitch and roll control. The inboard elevons 
are driven by a mechanical linkage commanded by the sum of pilot stick 
motion and the SAS (stability augmentation system) servo motion. The out- 
board elevon actuators are each slaved to the inboard elevon actuators by 
means of a mechanical linkage across the engine nacelles. 

The pitch-axis SAS used in the YF-12A is required full time and uses 
triple redundancy for reliability. A block diagram of this SAS is shown in 
Figure 2. The pitch SAS uses a pitch rate feedback for normal short-period 
stability augmentation and, at high altitude, also uses a lagged pitch rate for 
static stability augmentation. The gain of the SAS is scheduled with dynamic 
pressure as shown in Figure 2. 

The structural and aerodynamic data used in the study were provided for 
the longitudinal axis at the following flight conditions: 


Flight 

condition 

Mach 

Altitude, ft 

Weight, lb 

1 

0. 95 

25 000 

124 271 

1A 

0. 95 

25 000 

68 693 

2 

Supersonic 

Cruise 

90 703 

3 

0. 68 

22 000 

68 693 
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Figure 1. The YF-12A Aircraft 
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Figure 2. YF-12A Baseline SAS Block Diagram 













Analytical Representation 


The basic model used to analytically represent the YF-12A in the study 
•did not provide a rigorous description of the vehicle. However, the data 
used were readily available and were considered adequate for meeting the 
study objective. The basic model of the aircraft was comprised of a struc- 
tural model and an aerodynamic and gust model. 

Structural model . - The structural model included two-degree-of-free- 
dom rigid body dynamics and 16 symmetric structural modes. The model 
was capable of approximating the vehicle dynamics up to approximately 150 
rad/sec. Figure 3 shows typical mode shapes for the lower-frequency sym- 
metric structural modes. The model is further described in Appendix A. 

Aerodynamic and gust model . - The indicial aerodynamic formulation 
was used in this study to represent the unsteady aerodynamics in the time 
domain for use in the optimal control study. The aerodynamics are de- 
scribed in Appendix A. 

The atmospheric turbulence or gust model chosen for this study was the 
Press-Meadows model described in ref. 7. The effects of streamwise pene- 
tration into the gust were not included in the gust aerodynamic load formula- 
tion. This effect was considered minor due to the high velocity of the vehicle. 
The gust model is described in more detail in Appendix A. 

Data Modification . - The basic l8th-order model was modified by a 
change in variables to remove the zero roots. This was necessary to assure 
convergence of the computer program which computes the rms of a given air- 
craft parameter for a disturbance input. This change of variables is de- 
scribed in Appendix A. In the following text, the model resulting from the 
change of variables is called the "transformed 18th-order model. " 

The final step in the data modification consisted of generating a simpli- 
fied vehicle model for preliminary analysis studies. A simplified modelwas 
desired primarily to reduce the computation time required on the digital com- 
puter. The resulting simplified model was derived from the transformed 
1 8-degree-of-freedom model by neglecting structural modes 9 through 16. 
These higher-frequency modes had natural frequencies in excess of 100 rad/ 
sec and made a negligible contribution to the aircraft acceleration responses 
to a gust input. Although the simplified model was used in the preliminary 
analysis, the transformed 18th-order model was used to obtain a final per- 
formance and stability check. 

Hence, as a consequence of the data modification, three analytical 
representations of the aircraft were used in the study. These representa- 
tions are referred to as follows: 

« Basic 18th-order model 

• Transformed 18th-order model 

• Simplified model. 
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YF-12A response characteristics. - The basic response characteristics 
of the YF-12A were examined prior to the synthesis of LAMS control systems. 
The purpose of this task was threefold, namely: 

• To verify the transformed 18-degree-of-freedom aircraft 
model derived from the basic model 

• To verify the simplified model which was extracted from 
the transformed 18-degree-of-freedom model 

• To establish the nominal aircraft response characteristics. 

The transformed 18th-order model and the simplified model were veri- 
fied by comparison of eigenvalues and frequency responses, computed with 
these models, with similar data furnished for the basic model. Eigenvalues 
are compared in Table I for F/C 1, the high subsonic, heavyweight flight 
condition, with the baseline SAS loop opened. The data in Table I shows that 
the only differences in the eigenvalues of any significance are in the short- 
period mode. The differences shown were considered acceptable for the 
purposes of the study. 

Frequency responses computed with the transformed 18th-order model 
and the simplified model are compared with the basic 18th-order model re- 
sponses in Figure 4. These comparisons show a normal acceleration re- 
sponse at FS 248 to a combined inboard/ outboard elevon input disturbance 
for F/C 1. The response comparisons in Figure 4 show that above 1 rad / 
sec the responses are in good agreement. The differences at low frequency 
between the basic model and the transformed model are attributed to the 
absence of poles at the origin in the transformed model. The differences 
shown in Figure 4 were considered acceptable. The true test of the simpli- 
fied model was made when the systems designed using the simplified model 
were evaluated using the transformed model. These results are discussed 
later in the text, under SYSTEM PERFORMANCE. 

Power spectral density (PSD) plots of normal acceleration for a 1 -ft /sec 
gust input are shown in Figure 5 for the transformed 18th-order model. 

Plots are shown for acceleration responses at FS 248, 738, and 1236 for 
both F/C 1 and F/C 2. These data illustrate the performance characteristics 
of the nominal aircraft with the baseline SAS loop closed. The plots show the 
relative contribution of each of the aircraft modes to the acceleration levels 
of the selected aircraft stations. As these plots show, the rigid body and 
the first five or six modes were the primary contributors. 

A study of the plots in Figure 5 shows that total suppression of the 
structural mode contributions would reduce the rms acceleration by 40 per- 
cent at the pilot's station. On the other hand. Figure 5 shows only a small 
contribution to the rms acceleration by the structural modes at FS 738. 

At F/C 2, the responses in Figure 5 show the structural modes are the 
dominant contributors to the acceleration response. Total suppression of 
the structural modes would effect nearly a complete reduction in gust-induced 
accelerations . 
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TABLE I.- COMPUTED EIGENVALUES FOR F/C 1 



Basic 18th-order model 
without baseline SAS 

Transformed 18th-order 
model without SAS 

Simplified model 
without SAS 

Root 

Quantity 

Damping 

ratio 

Frequency, 
rad/ sec 

Quantity 

Damping 

ratio 

Frequency, 
rad/ sec 

Quantity 

Damping 

ratio 

Frequency, 
rad/ sec 

z, e 

2 

0 

0 

0 

Omitted 

— 

0 

Omitted 

— 

Rigid body 

2 

0. 37698 

1. 9927 

2 

0. 2878 

2. 2876 

2 

0. 3433 

1. 9929 

1st mode 

2 

0. 05232 

12. 4960 

2 

0. 05250 

12. 5070 

2 

0. 05570 

12.4076 

2nd mode 

2 

0. 05096 

16. 6216 

2 

0. 05087 

16.6515 

2 

0. 04 705 

16. 6126 

3rd mode 

2 

0. 04310 

29. 8090 

2 

0. 04306 

29. 8525 

2 

0. 04286 

29. 7795 

4th mode 

2 

0. 03946 

44. 6555 

2 

0. 03943 

44. 7127 

2 

0. 03941 

44. 5538 

5th mode 

2 

0. 04715 

49.6822 

2 

0. 04709 

49. 7628 

2 

0. 044 1 5 

49. 3676 

6th mode 

2 

0. 02823 

60. 7835 

2 

0. 02822 

60. 8390 

2 

0. 02808 

60. 8285 

7th mode 

2 

0. 02897 

73. 7769 

2 

0.02896 

73. 8450 

2 

0. 02831 

73. 8785 

8th mode 

2 

0. 03101 

86. 9578 

2 

0. 03100 

87. 0434 

2 

0. 03064 

87, 0373 

9th mode 

2 

0. 03077 

96. 8212 

2 

0. 03075 

96. 9162 

1 

Omitted 

— 

10th mode 

2 

0.03830 

105. 7804 

2 

0. 03827 

105. 9097 

I 

Omitted 

.... 

11th mode 

2 

0. 02991 

124. 7189 

2 

0. 02989 

124.8380 

EH 

Omitted 

.... 

12th mode 

2 

0. 03647 

131. 9372 

2 

0. 03644 

132. 0932 

0 

Omitted 

— 

13th mode 

2 

0. 03056 

146.7541 

2 

0. 03055 

146. 8977 

0 

Omitted 

— 

14 th mode 

2 

0. 03282 

150. 8634 

2 

0. 03280 

151. 0222 

0 

Omitted 

— 

15th mode 

2 

0. 02822 

161. 0100 

2 

0. 02821 

161. 1566 

0 

Omitted 

— 

16th mode 

2 

0. 03736 

156. 5301 

2 

0. 03734 

156. 7187 

0 

Omitted 

— 

Indie ial 
functions 

- 

Not listed 

— 

18 

— 

9. 77 

10 

— 

9. 77 

Series 

servo 

* 

Not listed 

— 

2 

0. 473 

133. 0 

2 

0.473 

133. 0 

Rate gyro 

- 

Not listed 

— 

2 

0. 602 

144. 0 

2 

0.602 

144. 0 

Inboard 

actuator 

- 

Not listed 

— 

2 

0.637 

39. 6 

2 

0. 637 

39. 6 

Outboard 

actuator 


Not listed 

— 

2 


33. 9 

2 

— 

33. 9 

Wind 

filter 

- 

Not listed 

— 

2 


0. 978 

2 

— 

0. 978 

Gust lag 

“ 

Not listed 


1 

— 

4. 88 

1 

— 

4. 88 
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Figure 4. Comparison of Frequency Responses Obtained With Basic 18th 
Order Model, Transformed 18th-Order Model, and Simplified 
Model, Without SAS 
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Figure 5A. Mean Square Acceleration at FSs 248, 738, and 1236 for 
a 1-ft/sec Turbulence Input -- Aircraft With SAS -- 
F/C 1 




Aircraft response with SAS 

• F/C 2 

• Transformed l8th-order model 
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Figure 5B. Mean Square Acceleration at FSs 248, 738, and 1236 for 
a 1-ft/sec Turbulence Input -- Aircraft With SAS -- 






Table II compares rms acceleration levels for a unit gust input at all 
four flight conditions. This table shows that F/Cs 1 and 1A have the largest 
rms acceleration levels. It is also of interest to note that for the three 
stations shown, the largest rms acceleration occurs at the aft end of the 
fuselage. 


TABLE II.- RMS ACCELERATION PER 1-FT/SEC GUST INPUT 
(WITH BASELINE SAS) 


Fuselage 

station 

F/C 1, 

high subsonic, 
heavyweight 
condition 

F/C 1A, 
high subsonic, 
lightweight 
condition 

F/C 2, 
supersonic 
cruise 
condition 

F/C 3, 

low subsonic, 
lightweight 
condition 

248 (pilot) 

0. 02118 g 

0. 02751 g 

0. OllOOg 

0. 01478 g 

738 

0. 01473 

0. 02173 

0. 00551 

0. 01492 

1236 (aft fuselage) 

0. 02683 

0. 04034 

0. 00996 

0. 01757 


LAMS System Synthesis 

Before discussing the detailed system design, it is worthwhile to con- 
sider some general aspects of the problem. These aspects include the de- 
sign criteria and design constraints, the use of force producers to achieve 
effective mode suppression, and the methods used to synthesize the systems. 

Design Criteria and Constraints 

A design goal in terms of structural mode suppression was established 
for the purposes of the study. This goal was to reduce the structural mode 
contribution to the pilot's rms acceleration by 50 percent for either a gust 
or control input disturbance. It was assumed that a 50-percent reduction 
would be sufficient for demonstration purposes. 

The systems were to be designed to meet this performance goal, but 
they also had to satisfy the following design constraints: 

• The LAMS system shall not degrade existing handling 
qualities provided by the baseline SAS. 

• The canards or mass-reaction device shall be restricted 
in size and force output so as not to cause a significant 
change in the YF-12 aerodynamics and to maintain internal 
loads within allowable limits. For the study, the design 
objective was to restrict the force output on the fuselage 
nose to less than 600 lb. 
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• The baseline SAS system operation shall not be altered and 
must remain engaged at all times. 

• The mechanical linkage connecting inboard and outboard 
actuators must be retained. 

• LAMS sensor locations must be physically practical. 

• The LAMS system shall meet stability margin requirements 
of a 6-dB gain margin and a 60-deg phase margin. 

LAMS Force Producers 

As stated above, the ground rules stipulated that the LAMS system could 
not alter the operation of the baseline SAS nor degrade the rigid-body hand- 
ling qualities. This requires that any coupling be minimized between the 
control of the aircraft's structural modes and the control of the rigid-body 
motion. The coupling between these two control functions may be minimized 
by either one or a combination of the following methods (assuming a prede- 
fined rigid-body SAS): 

• Suitably locating and sizing the LAMS force producers 

• Properly selecting the LAMS sensor locations 

• Properly designing the filtering in the LAMS system. 

The LAMS force producers must be selected before the sensor configu- 
ration and filtering can be established. To minimize coupling with rigid- 
body control, it is desirable to locate a LAMS force producer on the aircraft 
where structural mode deflections are large. This suggests that the force 
producers for LAMS should be located near the extremities of the aircraft, 
such as on the wing tips or on the nose. The only existing force producer on 
the YF-12A which qualifies (in the longitudinal axis) is the outboard elevon. 
But, on the YF-12A, the outboard elevon is slaved to the inboard elevon 
through a mechanical linkage. There is also no provision for any other in- 
puts (either mechanical or electrical) to the outboard actuator. Without in- 
dependent operation of the inboard and outboard elevons, the use of the force 
producer to separate the rigid-body control from the structural mode control 
would be ineffective. Hence, in this case, separation of rigid-body and 
structural mode control functions may be achieved only by a judicious loca- 
tion and blending of sensors and/or by filtering. Locating and blending of 
sensors is fairly effective, but this method is also limited by space available 
and environmental factors such as temperature variations. Minimizing 
coupling between rigid body and structural mode control by filtering is much 
less effective than either of the other two methods. Reliance on filtering to 
minimize the coupling generally results in a significant performance com- 
promise between the two functions. This is true especially when the rigid- 
body short-period frequency is close (less than a factor of 3) to the lowest 
structural mode frequency. It was concluded that if effective structural mode 
control was to be achieved on the YF-12A with the existing force producers. 
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the aircraft would have to be modified. This modification would be necessary 
to permit independent operation of the inboard and outboard elevons to en- 
hance the effectiveness of the outboard elevons for structural mode control. 

Because the outboard elevons are located near the wing tips, they should 
qualify as an effective LAMS force producer for controlling wing flexure. 
Control of the aircraft's fuselage modes, however, may be better achieved 
through the use of a force producer located on the nose of the aircraft. Ex- 
amination of the structural mode shapes (see Figure 3) shows the mode de- 
flections along the fuselage are largest at the nose and are relatively small 
at the tail (near the inboard elevons). It was concluded the inboard elevons 
would be relatively ineffective as a LAMS force producer because of potential 
coupling problems with the rigid-body control function. Because of these 
considerations, use of an auxiliary force producer located near the nose of 
the aircraft was considered in the study. Use of such a force producer would, 
of course, also require modification of the aircraft. 

Three primary factors must be considered before a force producer can 
be added. First, there must be space available within the aircraft for the 
force producer actuation. Second, the local structure must be capable of 
supporting the force producer. Third, the force producer cannot signifi- 
cantly alter the aerodynamics, especially with respect to the engine inlet 
dynamics. Small canard control surfaces or possibly a mass-reaction de- 
vice were considered as candidates for the auxiliary force producers. Their 
location on the aircraft would be as shown in Figure 6. 

Canard vanes were selected as one of the candidates because of promis- 
ing results obtained with them in the XB-70 program. Small canard vanes 
were installed on the XB-70 primarily for use as a means for exciting the 
structural modes. The flight testing performed on the XB-70 demonstrated 
their effectiveness for exciting the modes. If they were effective for exciting 
the modes, they should be effective for controlling the modes. Furthermore, 
because of their small size (4 sq ft) the XB-70 canard vanes had only a limi- 
ted effect on the rigid-body motion. One objective of the YF-12A study was 
to determine the size of canard vanes required for the YF-12A, keeping in 
mind the design constraints. The canard size required for mode control had 
to be evaluated against design constraints on allowable internal loads and on 
the aerodynamic effects on the engine inlets. 

The mass-reaction device was considered because of its potential ad- 
vantages over an aerodynamic control surface. One advantage is that it is 
not subject to the gross changes in aerodynamics encountered on vehicles 
like the YF-12A. This implies that gain scheduling of a system using mass 
reaction would be relatively simple. Because the mass is mounted internal 
to the aircraft, its addition would not alter the aerodynamic efficiency or 
cause potential engine inlet problems. 

Use of a mass-reaction device also has disadvantages. First, it appears 
to be somewhat self-defeating in that it adds weight to the aircraft. Second, 
it represents a greater design risk than does an aerodynamic surface. This 
risk stems from the fact that it has not been used in an airplane as a force 
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Figure 6. YF-12A Aircraft Showing Possible Locations 
of Auxiliary Force Producers 
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producer with feedback control. However, it was decided that the mass- 
reaction device should be evaluated to determine if the advantages outweigh 
the disadvantages. As with the canard, it was necessary to determine the 
size, displacement, and means for actuation of the mass. Finally, it was 
necessary to analyze these parameters for the YF-12A to determine if such 
a system would be reasonable for the space available and if the loads pro- 
duced would be within the allowed limits. Two principal assumptions were 
made concerning the mass-reaction force producer: 

1) It was assumed to be mounted at FS 248. A mass displacement 
of up to plus or minus 1 ft would be considered acceptable. 
Significantly larger displacements would have to be considered 
in more detail than intended for this study to determine the 
feasibility. At the bulkhead near FS 248, a space approximately 
4 ft high is available. Consequently, allowing for packaging of 
hardware, a plus or minus 1-ft displacement allowance appeared 
reasonable. 

2) The mass was assumed to be free to move in only a vertical 
plane and was to be positioned by an actuator (most likely 
hydraulic). In the absence of a command input to the actuator, 
the mass can be considered as rigidly attached to the airframe 
through the hydraulic actuator ram. 

To summarize, it was evident at the beginning of the study that some 
modifications of the aircraft would be required in order to demonstrate 
significant LAMS performance. The purpose of the modifications would be 
to provide some means of achieving structural mode control without de- 
grading rigid-body control. The modifications would consist of altering the 
outboard elevon servo to accept inputs from a LAMS system and/or adding 
an auxiliary force producer near the nose of the aircraft. The auxiliary 
force producer would be in the form of a small canard vane or a mass- 
reaction device. Combinations of all force producers were to be evaluated 
in the study to determine which combination would offer the greatest benefits. 

Synthesis Techniques 

A special frequency response technique was used as the basic design 
tool to configure LAMS controllers in the study. The technique was gener- 
ated during development work on the Lockheed SST ride quality problem 
(ref. 6). The primary motivation for its development at that time was to 
provide a means for synthesizing LAMS systems when vehicle data were 
provided only in a frequency response format. The general principles of 
the technique are described in Appendix B. The technique basically con- 
sists of a set of useful formulas and guidelines for specifying the LAMS 
sensor configuration and corresponding gains and filtering. The formulas 
specify requirements for control system gains and filtering in terms of a 
performance parameter and the imposed stability constraints. This tech- 
nique has proven to be a useful design tool. However, it also has some 
limitations. First, it becomes cumbersome to use in multiple-input (i. e. , 
multiple-force-producer) applications. Multiple-input systems designed 
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with this technique must be done through the usual procedure of designing 
one control loop at a time. Determining the allocation of control among the 
available force producers is an iterative trial and error procedure. The 
second limitation is that one must initially assume an allocation of mode 
suppression performance for each of the significant structural modes. 

Whether or not this performance allocation is optimum with regard to all 
factors, such as system complexity and handling quality constraints, cannot 
be determined except by trial and error. These are not serious problems 
for the experienced designer, but it would be desirable to establish a more 
^systematic procedure. 

Quadratic optimization from modern control theory (ref. 10) has been 
used in previous LAMS programs (ref. 4) as an alternate synthesis technique. 
The primary motivation for its use is the ease with which it can handle the 
multiple-input/multiple-output problem of which LAMS is typical. By proper 
selection of weighting factors in a performance cost function, the allocation 
of control among the force producers is automatically optimized. Further- 
more, quadratic optimization provides an effective means for determining 
the allocation of mode suppression performance among the structural modes. 

A limitation in applying quadratic optimization, however, is the fact that 
the feedbacks determined by the process are not specified in terms of prac- 
tical sensor outputs (e. g. , an accelerometer at a specified location). Rather, 
the feedbacks are specified directly in terms of the state variables. The 
transformation from these optimal feedbacks to practical sensor outputs re- 
mains a formidable task. 

One approach to "practicalizing" the optimal feedbacks is to use the fre- 
quency response technique. This is suggested by the fact that the two tech- 
niques appear to complement one another. Quadratic optimation provides an 
allocation of the control function among the given force producers but fails to 
define feedbacks in terms of practical sensor outputs. The frequency re- 
sponse technique, which lacks the capability for systematic allocation of the 
control function, offers the capability to define feedbacks in terms of practical 
sensors. A successful integration of these two techniques could yield a more 
effective synthesis procedure. An attempt was made in the YF-12A study to 
integrate these two synthesis techniques. The details of this effort, together 
with a description of the quadratic optimization technique, are presented in 
Appendix B. 

Quadratic optimization was successfully used in a previous program 
(ref. 4) to define a LAMS controller for the B-52. In that application, the 
optimal controller was transformed to a "practical" configuration by trial- 
and-error analysis on the analog computer. 

Design Procedure 

The steps of the design procedure used to define the mode suppression 
systems are described in the following paragraphs. 
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Three major steps were used. The first step was to assign a control 
function to each of the available force producers. The second step was to 
define the control system sensor/filter configuration. The third step was to 
evaluate system performance. This, of course, was an iterative process, 
which began with a set of aircraft data already in a suitable form for analysis 
and the available force producers specified. The three steps are discussed 
in more detail below. 

Step 1. Assign control function to each of the available force p roduc ers. - 
In this step, a control function was assigned to each force producer 7 . This 
assignment was based on the design objective, and the force producer's loca- 
tion and relative force output capability. The design objective was to mini- 
mize the contribution of structural modes to the acceleration response at 
selected fuselage stations. This was to be accomplished without degrading 
the rigid-body handling qualities. A force producer was considered ideal for 
mode suppression if it exhibited the following properties: 

• If it was located at or near a station where the acceleration 
response was to be reduced 

• If it possessed relatively large force coefficients for the 
significant structural modes 

• If it possessed relatively small rigid-body force coefficients. 

For example, a small canard vane located at the pilot's station would be a 
logical choice for mode suppression at the pilot's station. However, it would 
be relatively ineffective for suppressing the acceleration response on the 
wing tip unless the wing and fuselage flexure modes were strongly coupled. 

A wing tip-mounted force producer would be a better candidate for suppres- 
sing the structural mode acceleration response on the wing tip. 

Step 2. Define control system sensor/ filter configuration . - Definition 
of the sensor/ filter configuration was started by assuming an acceleration 
feedback from an accelerometer located near the force producer. This 
configuration generally guarantees that the sensed structural mode signals 
will be of the correct phase relative to one another to provide control of the 
modes through the given force producer (see ref. 3). After picking this 
initial sensor, the next task was to establish the extent to which the modes 
could be controlled, i. e. suppressed, taking into account other design objec- 
tives such as handling qualities. The rigid-body content of the accelerometer 
feedback signal generally degrades the handling qualities and, in particular, 
the short-period damping. As a consequence, the relationship between 
degraded handling qualities and mode suppression performance was evaluated. 
This was done by first defining a filter for the accelerometer feedback, using 
the frequency response technique (see Appendix B), and then by evaluating the 
closed-loop eigenvalues and PSD plots of acceleration for a gust input. 

The next task was to reduce the effect of acceleration feedback on rigid- 
body response by changing the sensor configuration. There were two basic 
alternatives to this step. One was to blend two or more accelerometers to 


21 



minimize the sensed rigid-body signal. The other alternative was to add a 
rate gyro feedback to restore the short-period damping. If the force pro- 
ducer is located aft of the center of gravity, blending of accelerometers is 
preferred. Mounting a single accelerometer aft of the center of gravity 
(necessary for bending-mode stability) degrades system performance because 
of the attendant nonminimum phase property. This degradation cannot be 
eliminated by adding a rate gyro feedback but can be minimized by blending 
accelerometers. Satisfactory blending of accelerometers can be achieved by 
placing an accelerometer on each wing, averaging their signals, and sub- 
tracting a signal from a fuselage-mounted accelerometer. Ideally, for this 
purpose, the fuselage-mounted accelerometer should be located the same dis- 
tance aft of the center of gravity as the wing-mounted accelerometer. 

If the force producer is located ahead of the center of gravity, then the 
second alternative is preferred, that is, to add a rate gyro feedback. Effec- 
tive cancellation of the rigid-body signal from a blend of accelerometers for 
a forward fuselage- located force producer generally requires at least three 
accelerometers. Blending of three accelerometers for this purpose is diffi- 
cult to achieve while at the same time maintaining favorable structural mode 
feedback signals. However, with a forward fuselage-located force producer 
the addition of a rate gyro feedback would restore the short-period character- 
istics. The problem with the rate gyro is to locate it so as to minimize un- 
favorable structural mode pickup by the gyro. Since the main function of the 
gyro is to augment the short period damping, its output could be filtered to 
further minimize its influence on the structural modes. The filtering required 
can be computed using the frequency response technique. 

The discussion so far has considered the use of only one force producer. 

If another force producer is available which is effective for controlling the 
rigid body, then this would offer yet another alternative for restoring the 
handling qualities without sacrificing the mode suppression performance. 

The next task in defining the sensor configuration was to analyze the 
effects of intermodal coupling and the attendant stability constraints. Inter- 
modal coupling is especially significant when one is attempting to suppress 
structural mode contributions to acceleration at one station by using a force 
producer located at some other station. For example, intermodal coupling 
would be a problem if one were attempting to suppress structural mode con- 
tributions to pilot's acceleration by using a wing tip-mounted force producer. 
Most of the structural mode contribution to pilot's acceleration arises from 
fuselage flexure. But a wing-mounted force producer has greater effective- 
ness for suppressing wing modes then for suppressing fuselage modes. The 
intermodal coupling can dilute the effectiveness of the wing-mounted force 
producer. The problem was analyzed in the study by using the frequency 
response technique of Appendix B. The technique was applied to specify 
allowable ratios of structural mode shape displacements for achieving mode 
suppression without sacrificing system stability due to intermodal coupling. 
These ratios served to characterize the required sensor configuration. 

Knowing the allowed ratios, combinations of sensors were sought which 
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satisfied this mode shape constraint on the sensor signal content. Having 
obtained a best "fit" sensor combination, a corresponding filter was defined. 

Step 3. Evaluate system performance . - In this step, the performance of 
the potential system configuration defined in Step 2 was evaluated in terms of 
the short-period frequency and damping, aircraft-system stability, and mode 
suppression performance. If the performance objectives were not met, the 
design procedure was repeated until the deficiencies were eliminated. 

Once an acceptable system was defined, additional performance param- 
eters were evaluated. This analysis included an evaluation of structural 
loads, steady-state stick gradients, and actuator characteristics. The 
effectiveness of each force producer as a means for introducing a calibrated 
test input for flight test evaluation was also evaluated. 

System Configurations 

The LAMS system design was started by examining the possible com- 
binations of force producers. The 11 possible combinations are shown in 
Table III. Five of the more promising combinations were selected for de- 
tailed study. They are: 

• System A - Inboard and outboard elevons (operating independently) 

• System B - Mass reaction only 

• System C - Mass reaction plus outboard elevons 

• System D - Canard only 

• System E - Canard plus outboard elevons. 

Of these five combinations, only one. System A, used the inboard ele- 
vons. This system was studied to establish the potential benefits of using 
only the force producers which presently exist on the aircraft. 


TABLE III. - POSSIBLE FORCE PRODUCER COMBINATIONS 


Additional force 
producer(s) 

Force Producer 

Canard 

Mass 

reaction 

Inboard 

elevon 

Outboard 

elevon 

None 

X 

X 

X 

X 

Outboard 

elevon 

X 

X 

X 

— 

Inboard 

elevon 

X 

X 

— 

— 

Inboard and 
outboard elevons 

X 

X 

— 

— 




The design of each of the five systems is described in the following para- 
graphs. For each system, a block diagram and a description of the function 
of each element in the system is presented. 


System A - Inboard elevon plus outboard elevon. - Figure 7 is an anal- 
ytical block diagram of System A. This system'was defined primarily to 
demonstrate the performance potential and problem areas associated with 
using only the existing force producers. One of the design constraints im- 
posed was that the mechanical linkage between the inboard elevon actuator 
and the outboard elevon series servo could not be altered. This mechanical 
linkage is the only means available for commanding outboard elevon deflec- 
tions. Hence, the outboard elevons are slaved to the inboard elevons. For 
LAMS control, however, it is desirable to have independent control of the 
outboard and inboard elevons to enhance control of the wing modes. Inde- 
pendent control of the two control surfaces would require modification of the 
outboard elevon series servo. The series servo would have to be changed to 
accept both the mechanical linkage input and an electrical input from the 
LAMS system. In addition, any LAMS input to the inboard elevon must be 
electrically subtracted from the outboard elevon. Figure 7 shows this can- 
cellation path. The dynamics shown on this path simulate the inboard actu- 
ator dynamics. A potential problem would exist with this cancelling of a 
mechanical signal path with an electrical equivalent. Any variations in the 
mechanical system (i. e. , the power actuator and mechanical linkage) due to 
nonlinearities (such as hysteresis) and temperature effects would reduce the 
effectiveness of the cancellation path. 

As Figure 7 indicates, the pitch SAS provides an input to only the inboard 
elevon actuator, while the LAMS system provides inputs to both the inboard 
and outboard elevon actuators. The primary function of the LAMS outboard 
elevon controller is to suppress the aircraft's wing bending modes. The 
LAMS inboard elevon controller functions to control the fuselage bending 
modes. 

The LAMS outboard elevon controller uses three accelerometers. Of 
these three, one is placed on each wing close to the outboard elevons (i. e. , 

WS 294). These two wing-mounted sensors are averaged together to provide 
a measurement of wing bending. A third accelerometer is located on the 
fuselage at station 954 (aft wheel well). The signal from this accelerometer 
is subtracted from the averaged signal of the wing-mounted sensors. The 
third accelerometer is required to minimize the net pickup of rigid-body 
motion and to assure proper phasing on the aeroelastic mode signals. It was 
desirable to minimize the rigid-body pickup to maintain acceptable handling 
qualities. Ideally, the rigid-body content of the wing-mounted sensor signals 
can be cancelled by locating the fuselage-mounted accelerometer at the same 
distance back from the center of gravity as are the wing-mounted accelerom- 
eters. Unfortunately, locating the fuselage- mounted accelerometer in this 
manner will produce unfavorable phasing on the pickup of higher-frequency 
aeroelastic modes. But by moving the accelerometer forward on the fuselage, 
an acceptable compromise can be achieved between rigid-body pickup and 
phasing of the aeroelastic mode signals. The phasing on the aeroelastic mode 
signals was most critical for the first two aeroelastic modes. A strong inter- 
modal coupling exists between these two modes because of their close proximity 
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Figure 7. LAMS Inboard /Outboard Elevon Controller 


in frequency. The first mode is primarily a fuselage-bending mode while the 
second mode is primarily a wing-bending mode. The intermodal coupling 
aggravated the design of a LAMS system because, as control was applied to 
suppress the wing mode via the outboard elevons, the response of the fuselage 
mode was made worse. The problem was solved using the frequency re- 
sponse technique (Appendix B) to formulate a design constraint on the required 
modal pickup and the corresponding filter. The result was the outboard ele- 
von system shown in Figure 7. The filter which was computed provides 
approximately a 90- deg phase lag over the frequency range of the significant 
aeroelastic modes. The gain value of 0. 01 rad/g was computed to maximize 
modal suppression performance compatible with stability constraints. 

The inboard elevon system uses a blend of four accelerometers mounted 
along the fuselage. This sensor combination is used to eliminate rigid-body 
pickup and to provide the required phasing on the structural modes. Two 
accelerometers are required at FS 1236, since it was assumed a single 
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accelerometer could not be physically mounted on the centerline. The aver- 
aged output of these two aft-mounted accelerometers is blended with an accel- 
erometer mounted just aft of the cockpit at station 388 and an accelerometer 
mounted in the wheel well at station 954. The gains and shaping shown in 
Figure 7 for the inboard elevon system provide the required compensation f or 
augmentation of the fuselage bending modes, especially modes 1 and 3. 

The synthesis of the inboard elevon system presented the most difficult 
task of all the control loops considered. The basic difficulty with using this 
force producer for mode suppression is that the mode deflections are rela- 
tively small in the vicinity of the inboard elevons. Hence, the coupling 
problem with the rigid body was more serious than with the other force pro- 
ducers. Use of an accelerometer near the inboard elevons was ruled out be- 
cause of the potential degradation in rigid-body dynamics. Various combina- 
tions of rate gyros and accelerometers were considered, with the objective 
of providing favorable pickup on all structural modes while minimizing rigid - 
body pickup. The end result was the sensor configuration shown in Figure 7. 
This resulted in favorable pickup of the structural modes while eliminating 
the rigid-body content. 

Although this sensor configuration was considered acceptable for the 
purposes of the study, it possesses a potentially undesirable feature. This 
feature has to do with the blending of three (in effect) sensors distributed 
along the fuselage. Such a configuration is generally sensitive to changes 
in mode shapes, since proper blending for some of the modes will depend on 
taking the difference between two large numbers. 

System B - Mass reaction only . - Figure 8 is a block diagram of the 
mass-reaction control system. Since this force producer, like the canard, 
does not presently exist on the aircraft, the actuator dynamics were not 
specified. A brief analysis of the mass-reaction dynamics was made prior 
to synthesizing the mass-reaction system. The objective of this analysis 
was to establish the desired actuator frequency response characteristic for 
the mass-reaction system. The actuator response characteristic was deter- 
mined by examining the aircraft's acceleration response at the mass location 
for a mass -displacement input. This plot, shown in Figure 9, indicated that 
a low-frequency, second-order actuator characteristic would be desirable. 
Such an actuator would provide the following benefits: 

• A mass -position feedback 

• A 180-deg phase lag over the frequency range of the 
structural modes 

• A high-frequency gain attenuation. 

The mass-position feedback is required to keep the mass centered at 
least on a low-frequency basis. Without this feedback, the mass would be 
free to drift and eventually hit its displacement limits. A second-order 
actuator with a natural frequency of 2 rad/sec and a damping ratio of 0. 5 
was initially selected. An alternative to using this actuator would be to use 
a low-frequency, first-order actuator. The 180-deg phase lag required 
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would have to be obtained through a phase reversal in the control system. 
However, this approach was discarded in the study because the high-frequency 
attenuation characteristic was considered to be inadequate. 

The sensor configuration shown in Figure 8 consists of four accelerom- 
eters located along the fuselage. One accelerometer is located in the nose at 
FS 160. Another is located at FS 97 5, which is just aft of the baseline SAS 
rate gyro. The remaining two accelerometers are placed at FS 1212 on 
either side of the centerline. As was the case with the outboard elevon con- 
troller, it was assumed there would be no space available on the centerline 
for a sensor at this location. Consequently, the two accelerometers, placed 
on either side of the centerline, are averaged together. The four sensors 
were blended together to eliminate all rigid-body pickup while maintaining 
favorable structural mode pickup. As with the inboard elevon sensor con- 
figuration, blending of these accelerometers tends to result in a system which 
is not tolerant to variations in aircraft dynamics. This intolerance is dem- 
onstrated in the performance results. The gain and shaping shown in Figure 8 
was defined using the frequency response technique. 

System C - Mass reaction plus outboard elevon . - A block diagram of 
System C is shown in Figure 10. The mass-reaction controller shown in 
this figure is identical to the one described for System B. The outboard 
elevon controller shown is identical to the one used in System A except for 
the following: 

• Gain on the net sensor signal is doubled from 0. 01 rad/g 
to 0. 02 rad/g. 

• The blend of sensors shown calls for using an acceler- 
ometer mounted at FS 795. In System A, this sensor 
was mounted at station 955. This change was made to 
obtain commonality with the sensors used for the mass- 
reaction system. 

System C contains a rather complex combination of sensors. Refinement of 
this system would be directed toward simplifying this sensor complement. 

It is anticipated the two sensors at station 1212 could be eliminated. 

System D - Canard only. - An analytical block diagram of System D is 
shown in Figure 11. It is the simplest of all five candidate systems. The 
system requires use of a conventional (i. e. , no temperature compensation) 
accelerometer located near the canard at FS 248, and a signal from the base- 
line SAS rate gyro. The baseline SAS rate gyro is located on the fuselage 
(FS 765) approximately where the leading edge of the wing joins the fuselage. 
The only filtering required by the system is provided by the 0. 25-sec first- 
order lag specified for the actuator dynamics. 

With the canard located at the pilot's station, synthesis of a mode sup- 
pression system is relatively straightforward. The two primary design 
constraints in the synthesis are retention of handling qualities and system 
stability. The effects of these constraints are minimized by using an 
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Figure 10. LAMS Outboard Elevon-Plus-Mass Reaction Controller 



Figure 11. LAMS Canard-Only Controller 
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accelerometer placed at the same location as the canard. Use of an accel- 
erometer placed in this location offers the following advantages: 

• The structural mode pickup is large relative to other vehicle 
stations. This property minimizes any potential coupling 
problem with the rigid-body mode. 

• Placement of an accelerometer at the same location as the 
force producer tends to assure proper phasing for stability 
on all the structural modes. (This concept of locating the 
sensor at the same location as the force producers is de- 
scribed in ref. 3) 

• Phasing on the rigid-body acceleration pickup is such to 
provide some rigid-body load alleviation. 

The acceleration feedback to the canard results in commanding a positive 
canard deflection (i. e. , trailing-edge down) for a positive normal accelera- 
tion (i. e. , downward). A positive canard deflection results in a nose-up 
motion to reduce the angle of attack, thereby providing load alleviation. The 
price paid, however, is reduced short-period damping. The reduction in 
short-period damping is compensated by adding a pitch rate feedback to the 
canard. 

System E - Canard plus outboard elevons . - Figure 12 shows a block 
diagram of System E. In this system the outboard elevon loop configured 


Actuator 



Figure 12. LAMS Outboard Elevon-Plus-Canard Controller 
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for System A was added to the canard controller of System D. This was 
done to provide additional mode suppression of the wing aeroelastic modes. 

System Mechanization Considerations 

Mechanizations for the five candidate systems were examined in the 
study to the extent necessary to uncover problem areas affecting feasibility 
and to determine required aircraft modifications. The study indicated no 
insurmountable problems with respect to mechanization for any of the five 
candidate systems. All aircraft modifications required to incorporate the 
auxiliary force producers or to modify the existing force producers were 
reasonable. Installation of accelerometers in the wings or in the aft body 
would be constrained by the severe temperature environment, however. 
Either special high-temperature sensors would have to be designed, or 
cooling would have to be provided for conventional sensors or else the 
flight envelope would have to be restricted. 


SYSTEM PERFORMANCE 


Performance results are presented for each of the five candidate LAMS 
control systems. These results consist of the following vehicle/system 
properties : 

• Eigenvalues 

• Stability margins 

e Mode suppression performance at FS 24 8 (pilot's station), 

FS 738, and FS 1236 (aft fuselage) for a turbulence input 
disturbance 

• Mode suppression performance at the above three fuselage 
stations for input disturbances applied via a force producer. 

The force producer input disturbance was generated by 
assuming a white noise input into the gust model which was 
modified by setting L/V o equal to 20. [See eq. (A34) of 

Appendix A, Basic Vehicle Data]. The gust model was 
selected as a convenient means for attenuating the high- 
frequency content of the white noise input. The output of 
the modified gust model was summed into the control 
surface actuator being used to excite the aircraft. For 
the cases where the canard was used to excite the air- 
craft, actuator dynamics were ignored. 


Handling quality performance and the effects of nonlinearities are 
described collectively for the five systems at the end of the section. 
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System A - Inboard/Outboard Elevon System 


Table IV shows a comparison of the open- and closed-loop roots obtained 
with the inboard/outboard elevon system. Examination of these roots shows 
that the rigid-body damping was reduced and the frequency increased. The 
resultant frequency and damping were considered acceptable, however. The 
excessively high damping obtained on the first structural mode indicates the 
mode is being overcontrolled. These results indicate the gain on the inboard 
elevon controller should have been lower. If the gain were halved there 
would be less of a change in the short-period dynamics and probably no sig- 
nificant change in the mode suppression performance for the first mode. The 
mode suppression performance of the higher-frequency modes would have 
been reduced, but this could probably be restored via a change in the out- 
board elevon system. 


TABLE IV. - COMPARISON OF OPEN- AND 
CLOSED -LOOP POLES WITH SYSTEM A - 
F/C 1 - SIMPLIFIED MODEL 



Baseline SAS only 

With LAMS 

Mode 

Damping 

ratio 

Frequency, 
rad/ sec 

Damping 

ratio 

Frequency, 
rad/ sec 

Short-period 

0. 610 

3. 70 

0. 443 

5. 76 

1st mode 

0. 068 

12. 7 

0. 969 

11. 64 

2nd mode 

0. 056 

16. 2 

0. 109 

16. 2 

3rd mode 

0. 076 

29. 1 

0. Ill 

28. 0 

4 th mode 

0. 042 

44. 1 

0. 097 

44. 1 

5th mode 

0. 043 

49.4 

0. 058 

49. 1 

6 th mode 

0. 029 

60. 8 

0. 038 

60. 8 

7 th mode 

0. 032 

74. 1 

0. 044 

74. 3 

8th mode 

0. 029 

86. 8 

0. 074 

86. 1 


The stability results for the inboard/outboard elevon system are shown 
in Table V. Stability margins shown for F/C 1 tend to support the above ob- 
servation concerning the inboard elevon control gain. Table V indicates only 
a 4-dB margin exists at F / C 1. Reducing the inboard elevon control gain 
level would, of course, enhance this margin. The control gains were not 
changed for the computation of stability at F/C 1A. As the table indicates, 
the inboard elevon system causes a stability problem with the second mode at 
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TABLE V. - STABILITY MARGINS - SYSTEM A - 
TRANSFORMED 18TH-ORDER MODEL 


Flight 

Inboard elevon loop closed 

Outboard elevon loop closed 

condition 

Gain margin, dB 

Phase margin, deg 

Gain margin, dB 

Phase margin, deg 

1 

— 

180 

4 

63 

1A 

Unstabl 

e 2nd mode 

10 

63 

2 

Unstable higher- 

14 

67 


frequency mode 




this flight condition. However, closure of both the inboard and outboard con- 
trol loops restored the stability of this mode. The stability problem with the 
second mode is most likely due to intermodal coupling between the first and 
second modes. It is expected that a reduced gain for the inboard elevon loop 
would solve the problem. Stability at F/C 2 was computed with the inboard 
elevon gain increased by a factor of 10 and the outboard elevon gain increased 
by a factor of 6. These gain adjustments were made to compensate for the 
decrease in surface effectiveness in going from F/C 1 to F/C 2. As is evi- 
dent from the table, the inboard elevon controller causes an instability in a 
higher-frequency structural mode. This could be easily corrected by a gain 
reduction or by the addition of some high-frequency filtering. 

Table VI shows mode suppression performance obtained with the inboard 
plus outboard elevon system at two flight conditions. These results show a 
reduction in rms acceleration of approximately 30 percent at the subsonic 


TABLE VI. - MODE SUPPRESSION PERFORMANCE FOR A 
1 -FT/SEC TURBULENCE INPUT DISTURBANCE WITH 
SYSTEM A - TRANSFORMED 18TH-ORDER MODEL 


Fuselage 

station 

F/C 1, 

high subsonic, 
heavyweight condition 3 

F/C 2, 
supersonic 
cruise condition* 3 

SAS only 

LAMS 

SAS onlv 

LAMS 

248 (pilot) 

0. 0212 g 

0. 0159g 

0. 0110 g 

0. 0051 g 

738 

0. 0147 

0. 0142 

0. 0055 

0. 004 0 

1236 (aft fuselage) 

0. 0268 

0.0189 

0. 010 

0. 0050 


Simplified model 
^Transformed 1 8th-order model 
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condition and 50 percent at the supersonic condition. It should be noted that 
these rms values are misleading. The contribution of the rigid-body mode, 
which was intentionally not reduced, masks the true mode suppression per- 
formance. The corresponding PSD plots of mean square acceleration at the 
subsonic condition are shown in Figures 13, 14, and 15. The mode suppres- 
sion performance obtained is more evident from these plots. 

Tables VII and VIII show rms accelerations obtained for force producer 
input disturbances at the subsonic and supersonic conditions, respectively. 
Performance at the subsonic condition shows rms acceleration reductions up 
to 70 percent for an inboard/outboard elevon input. Performance obtained at 
the supersonic condition (see Table VIII) shows reductions in rms accelera- 
tions on the order of 50 percent for a canard input and 7 5 percent for an in- 
board/outboard elevon input. 


TABLE VII. - RMS ACCELERATION RESPONSE FOR A 
1-RAD FORCE PRODUCER INPUT DISTURBANCE - 
SYSTEM A - F/C 1 - SIMPLIFIED MODEL 


Fuselage 

station 

Input to 

canards 

Input to inboard/ 
outboard elevons 

SAS only 

LAMS 

SAS only 

LAMS 

248 

6. 52 g 

6. 16 g 

2. 16 g 

0. 82 5 g 

738 

1. 79 

1. 75 

0. 576 

0. 337 

1236 

2. 91 

2.23 

1. 96 

0. 586 


TABLE VIII. - RMS ACCELERATION RESPONSE FOR A 
1-RAD FORCE PRODUCER INPUT DISTURBANCE - 
SYSTEM A - F/C 2 - TRANSFORMED 18TH- 
ORDER MODEL 


Fuselage 

station 

Input to canards 

Input to inboard/ 
outboard elevons 

SAS only 

LAMS 

SAS only 

LAMS 

248 

0. 62 1 g 

0. 493g 

0. 52 6 g 

0. 134g 

738 

0. 220 

0. 167 

0. 222 

0. 076 

1236 

0. 300 

0. 149 

0. 376 

0. 106 
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Figure 13. Mean Square Acceleration at FS 248 for a 1-ft/sec Turbulence 
Input With and Without LAMS System A -- F/C 1 
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Figure 14. Mean Square Acceleration at FS 738 for a 1-ft/sec Turbulence 
Input With and Without LAMS System A -- F/C 1 
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Figure 15. Mean Square Acceleration at FS 1236 for a 1-ft/sec Turbulence 
Input With and Without LAMS System A -- F/C 1 





Table IX shows rms elevon displacements and rates obtained for each of 
the disturbances inputs. Surface displacements obtained for a turbulence in- 
put indicate that only 0. 3 deg of inboard elevon and 0. 11 deg of outboard ele- 
von rms displacement would be required for a 3a (12-ft/sec) gust. These 
values represent only the LAMS system surface displacements. Actual sur- 
face displacements would also include the contribution due to the baseline 
SAS command. Hence, it cannot be readily assessed from these results 
whether or not hysteresis in the actuator dynamics will seriously degrade 
the LAMS performance with this system. But since these rms values are 
indeed small, it is expected that hysteresis may pose a significant problem 
for turbulence inputs. 


TABLE IX. - RMS SURFACE MOTION FOR EACH DISTURBANCE 
INPUT - SYSTEM A - F/C 1 - SIMPLIFIED MODEL 


Disturbance 

input 

RMS inboard elevon 

RMS outboard elevon 

Displacement, 

rad 

Rate, 
rad/ sec 

Displacement, 

rad 

Rate 
rad/ sec 

1-ft/sec turbulence 
input 

0. 473 x 10* 3 

Not 

available 

0. 973 x 10' 4 

0. 00243 

1 rad input to 
canards 

0. 085 

Not 

available 

0. 00912 

0. 349 

1 rad input to in- 
board/outboard 
elevens 

0. 0275 

1.0 

0. 0078 

0.218 


The rms surface displacements and rates obtained for the input distur- 
bances applied to the force producers serve only to size the allowable level 
of input disturbance. The results obtained indicate no problems exist for 
reasonable rms acceleration levels. 


System B - Mass Reaction Only 

Table X shows the eigenvalues obtained with the mass-reaction system. 
These roots indicate the system augmented the lower-frequency fuselage 
modes (modes 1 and 3) without affecting the primary wing mode (mode 2). 
This illustrates the ineffectiveness of a nose-mounted force producer for 
suppressing wing modes. 

Table XI shows the stability margins computed at each flight condition 
with the transformed 18th-order model. The gain of the mass-reaction 
system was held constant for all three conditions. At F/C 1A the third mode 
was unstable. This situation was caused by a sign reveral on the force co- 
efficient for this mode. This suggests the mass should be moved forward to 
improve its tolerance to variations. A second alternative would be to change 
the sensor configuration. 
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TABLE X. - COMPARISON OF OPEN- AND CLOSED - 
LOOP POLES WITH SYSTEM B-F/Cl- 
SIMPLIFIED MODEL 



Baseline SAS only 

With LAMS 

Mode 

Damping 

ratio 

Frequency, 
rad/ sec 

Damping 

ratio 

Frequency, 
rad/ sec 

Short-period 

0. 61 

3.70 

0. 74 

2. 94 

1st mode 

0. 068 

12. 7 

0. 46 

10. 2 

2nd mode 

0. 056 

16. 2 

0. 056 

16. 3 

3rd mode 

0. 076 

29. 1 

0. 157 

28. 0 

4th mode 

0. 042 

44. 1 

0. 045 

44. 0 

5th mode 

0. 043 

49. 4 

0. 043 

49. 4 

6th mode 

0. 029 

60. 8 

0. 058 

60. 2 

7th mode 

0. 032 

74. 1 

0. 04 2 

73. 7 

8th mode 

0. 029 

86. 8 

0. 034 

86. 7 

Actuator 

0. 5 

4. 0 

0. 277 

4.47 


TABLE XI. - STABILITY MARGINS - 
SYSTEM B - TRANSFORMED 
18 TH-ORDER MODEL 


Flight 

Condition 

Gain margin, dB 

Phase margin, deg 

1 

9 

65 

1A 

Unstable 3rd mode 

Unstable 3rd mode 

2 

21 

55 





Table XII shows mode suppression performance for both turbulence and 
force producer inputs. This performance was computed at the high subsonic, 
heavyweight condition. Performance was not computed at the supersonic 
condition. It was felt that the canard performance obtained at the supersonic 
condition would suffice to show the feasibility of an auxiliary force producer. 
Figures 16, 17, and 18 show the PSD plots for a turbulence input. These 
plots show that the mass-reaction system is fairly effective for suppressing 
the lower-frequency modes. However, the system provides little benefit at 
frequencies above 6 Hz for a gust input. The mass- reaction system provides 
little improvement for the wing modes, especially at station 1236. 


TABLE XII. - RMS ACCELERATION RESPONSE TO TURBULENCE 
AND FORCE PRODUCER INPUTS - SYSTEM B - F/C 1 - 

SIMPLIFIED MODEL 


Fuselage 

station 

1-ft/sec turbulence 
input 

1-in. input to 
mass reaction 

1-rad input to inboard/ 
outboard elevons 

SAS only 

LAMS 

SAS only 

LAMS 

SAS only 

LAMS 

248 (pilot) 

0. 0211 g 

0. 0164 g 

1. 22 x 10- 4 g 

0. 911 x 10‘ 4 g 

2. 16 g 

0. 897 g 

738 

0. 0148 

0. 0146 

0. 334 x 10' 4 

0. 220 x 10' 4 

0. 576 

0. 397 

1236 (aft fuselage) 

0. 0253 

0. 0232 

0. 550 x 10‘ 4 

0. 364 x 10" 4 

1. 96 

1. 20 


Table XIII shows the rms values of mass displacement, rate, accelera- 
tion, and force output computed for a turbulence input. These results indi- 
cate that a 100-lb weight will require approximately a 1-ft displacement to 
produce the performance described in the previous paragraphs. If necessary, 
this displacement can be reduced by a proportionate increase in the weight of 
the mass. 


TABLE XIH. - MASS VARIABLES FOR A TURBULENCE 
INPUT DISTURBANCE - SYSTEM B - F/C 1 - 
SIMPLIFIED MODEL 


RMS mass 
variables 3 

Unit turbulence 
input, 1 ft/sec 

3 C7 turbulence input, 
12 ft/sec 

Displacement 

1. 1 in. 

13. 2 in. 

Rate 

7. 44 in. / sec 

89. 3 in. / sec 

Acceleration 

89. 0 in. 1 sec^ 

107 0 in. / sec^ 

Force output a 

23. 0 lb 

277 lb 


a Assumes 100-lb weight for the mass. 
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Frequency, Hz 

Figure 16. Mean Square Acceleration at FS 248 for a 1-ft/sec Turbulence 
Input With and Without LAMS System B -- F/C 1 
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Figure 17. Mean Square Acceleration at FS 738 for a 1-ft/sec Turbulence 
Input With and Without LAMS System B -- F/C 1 
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Figure 18. Mean Square Acceleration at FS 1236 for a 1-ft/sec Turbulence 
Input With and Without LAMS System B -- F/C 1 



Table XIII shows that the force output on the nose was 277 lb for a 3b 
gust input. This is well within the assumed force constraint of 600 lb. 

System C - Mass Reaction Plus Outboard Elevon 

The open- and closed-loop roots are compared in Table XIV for the 
mass-reaction-plus-outboard elevon system. These roots show significant 
augmentation of the first three structural modes. These results, when com- 
pared with the results obtained with the mass- reaction system only, indicate 
the effectiveness of the outboard elevon system to control the wing modes. 


TABLE XIV. - COMPARISON OF OPEN- AND CLOSED-LOOP 
POLES WITH SYSTEM C - SUBSONIC, HEAVYWEIGHT 
CONDITION - SIMPLIFIED MODEL 



Baseline SAS only 

With LAMS 

Mode 

Damping 

ratio 

Frequency, 
rad/ sec 

Damping 

ratio 

Frequency, 
rad / sec 

Short-period 

0. 610 

3. 70 

0. 743 

2. 94 

1st mode 

0. 068 

12. 7 

0. 418 

10. 6 

2nd mode 

0. 056 

16. 2 

0. 217 

16.4 

3rd mode 

0. 076 

29. 1 

0. 203 

27. 8 

4th mode 

0. 042 

44. 1 

0. 083 

45. 0 

5th mode 

0. 043 

49.4 

0. 096 

47. 7 

6th mode 

0. 029 

60. 8 

0. 059 

60. 2 

7th mode 

0. 032 

74. 1 

0. 04 2 

73. 7 

8th mode 

0. 029 

86. 8 

0. 034 

86. 7 

Actuator 

0. 5 

4. 0 

0. 267 

4.46 


Table XV shows the stability margins obtained with this control system 
for three flight conditions. The results show a stability problem exists with 
the third mode at F/C 1A. This instability results from closure of the mass- 
reaction control loop. This same problem was evident with the mass-reaction 
system only. It should be noted, however, that with both loops closed the 
system was stable. The outboard elevon loop compensated for the destabili- 
zing effect of the mass-reaction system. The control gains were held con- 
stant between conditions 1 and 2, but the outboard elevon gain was increased 
by a factor of 6 at F/C 2. 


44 




TABLE XV. - STABILITY MARGINS - SYSTEM C - 
TRANSFORMED 18TH-ORDER MODEL 


Flight 

Mass reaction loop closed 

Outboard elevon loop closed 

condition 

Gain margin, dB 

Phase margin, deg 

Gain margin, dB 

Phase margin, deg 

i 

— 

88 

9 

64 

1A 

Unsta 

3le 3rd mode 

3 

87 

2 

— 

102 

10 

57 


Table XVI shows the rms acceleration performance obtained for both 
gust and force producer inputs. PSD plots for a turbulence input disturbance 
are shown in Figures 19. 20, and 21. As expected, this rms acceleration 
performance is similar to the mass -reaction- only performance. The major 
differences to be noted are for the improved rms acceleration performance 
obtained at station 1236. This additional performance is attributed to sup- 
pression of the wing modes by the outboard elevon system. Performance of 
this system was computed for only the high subsonic, heavyweight condition. 


TABLE XVI. - RMS ACCELERATION RESPONSE TO 
TURBULENCE AND FORCE PRODUCER INPUTS - 
SYSTEM C - F/C 1 - SIMPLIFIED MODEL 


Fuselage 

station 

1 ft/sec turbulence 
input 

1-in. input to 
mass reaction 

1-rad input to inboard/ 
outboard elevons 

SAS only 

LAMS 

SAS only 

LAMS 

SAS only 

LAMS 

248 

0. 0211 g 

0. 0162 g 

1. 22 x 10' 4 g 

0. 902 x 10 _4 g 

2. 16 g 

0. 729 g 

738 

0. 0148 

0. 0147 

0. 334 x 10‘ 4 

0. 217 x 10" 4 

0. 576 

0. 338 

1236 

0. 0253 

0. 0222 

0. 550 x 10‘ 4 

0. 361 x 10' 4 

1. 96 

0. 833 


Table XVII shows rms values for the mass system variables obtained 
for a turbulence input. These results show the mass variables were in- 
creased by approximately 10 percent with the addition of the outboard elevon. 
The rms mass displacement of 14. 6 in. obtained with a 3a gust is larger than 
desired, although the force level remains well below the allowable force of 
600 lb. As noted previously, the displacement of the mass can be reduced by 
increasing the weight of the mass. 
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Figure 19. Mean Square Acceleration at FS 248 for a 1-ft/sec Turbulence 
Input With and Without LAMS System C -- F/C 1 
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Figure 20. Mean Square Acceleration at FS 738 for a 1-ft/sec Turbulence 
Input With and Without LAMS System C -- F/C 1 
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Frequency, Hz 

Figure 21. Mean Square Acceleration at FS 1236 for a 1-ft/sec Turbulence 
Input With and Without LAMS System C -- F/C 1 





TABLE XVII. - MASS VARIABLES FOR A 
TURBULENCE INPUT DISTURBANCE - 
SYSTEM C - F/C 1 - SIMPLIFIED 
MODEL 


RMS Mass 
variables 4 

Unit turbulence 
incut, 1 ft/sec 

3ct turbulence input, 
12 ft /sec 

Displacement, in. 

1. 22 

14. 6 

Rate, in. / sec 

9. 13 

110 

Acceleration, 
in. / sec^ 

102 

1220 

Force output, lb 

26. 5 

318 


a Assumes 100-lb weight for the mass. 


System D - Canard-Only Controller 

Table XVIII shows the eigenvalues obtained with the canard-only con- 
troller at F/C 1. Not too surprising is the fact that the damping of the 
second mode, a wing mode, is not significantly affected. As will be shown, 
however, its contribution to pilot's acceleration is significantly reduced. 

Table XIX shows the stability margins of this system at the two high 
subsonic, heavyweight flight conditions (1 and 1A) and at the supersonic 
cruise flight condition (F/C 2). These stability margins were extracted 
from open-loop frequency responses computed using the transformed 18th- 
order model. At F/C 2, the control gain in the canard system was increased 
by a factor of 10 over the one used at F/C 1 (or F/C 1A). 

Table XX shows rms acceleration performance obtained at the three 
fuselage stations for 1-ft/sec turbulence input. These results were com- 
puted using the transformed l8th-order vehicle model. The canard-only 
system is the only one for which performance was computed at the high 
subsonic, lightweight condition. It is also the only configuration for which 
performance was computed at all conditions with the transformed 18th-order 
model. The performance shown in Table XX shows significant reductions 
in rms acceleration at the pilot's station for all three conditions. A 40-per- 
cent reduction was obtained at the subsonic conditions and a 62-percent 
reduction was obtained at the supersonic condition. This performance is 
especially significant when the relatively large contribution of the rigid-body 
mode to rms acceleration is considered. As the subsonic, heavyweight con- 
dition the rigid-body mode contributes approximately 60 percent to rms 


49 




TABLE XVIII. - COMPARISON OF OPEN- AND 
CLOSED -LOOP POLES WITH SYSTEM D - 
F/C 1 - SIMPLIFIED MODEL 


Mode 

Baseline SAS only 

With LAMS 

Damping 

ratio 

Frequency, 
rad/ sec 

Damping 

ratio 

Frequency, 
rad /sec 

Short-period 

0. 610 

3. 70 

m 

4. 24 

1st mode 

0. 068 

12. 7 


9. 89 

2nd mode 

0. 056 

16. 2 

m/BM 

16. 3 

3rd mode 

0. 076 

29. 1 

0. 141 

25. 9 

4th mode 

0. 042 

44. 1 

0. 059 

43. 2 

5th mode 

0. 043 

49. 4 

0. 044 

49. 4 

6th mode 

0. 029 

60. 8 

0. 067 

59. 8 

7th mode 

0. 032 

74. 1 

0.038 

73. 9 

8th mode 

0. 029 

86. 8 

0. 032 

86. 7 


TABLE XIX. - STABILITY MARGINS - 
SYSTEM D - TRANSFORMED 18TH- 
ORDER MODEL 


Flight 

condition 

Gain 

margin, 

db 

Phase 

margin, 

deg 

1 

10 

90 

1A 

— 

88 

2 

10 

78 


acceleration at the pilot's station. Figures 22 and 23 show the PSD plots of 
pilot's acceleration at the high subsonic, heavyweight and supersonic cruise 
conditions. 
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Figure 22. Mean Square Acceleration at FS 248 for a 1-ft/sec Turbulence 
Input With and Without LAMS System D -- F/C 1 
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Figure 23. Mean Square Acceleration at FS 248 for a 1-ft/sec Turbulence 
Input With and Without LAMS System D -- F/C 2 





The data in Table XX also show the relative ineffectiveness of the canard 
in providing mode suppression at the aft end of the fuselage. The inability to 
provide mode suppression at station 1236 is attributed to two factors. First, 
it was a design constraint not to provide rigid-body load alleviation with the 
canard. This constraint was imposed to minimize canard size requirement. 
Second, because the canard is located on the nose, it is ineffective for con- 
trolling wing modes, such as the second structural mode. Figures 24 through 
27 show the performance at both stations 738 and 1236 for the high subsonic, 
heavyweight and the supersonic conditions. 


TABLE XX. - MODE SUPPRESSION PERFORMANCE FOR 
A 1 -FT/SEC TURBULENCE INPUT DISTURBANCE - 
SYSTEM D - TRANSFORMED 1 8TH- ORDER MODEL 


Fuselage 

station 

F/C 1, 

high subsonic, 
heavyweight 
condition 

F/C 1A, 
high subsonic, 
lightweight 
condition 

F/C 2 , 
supersonic 
cruise 
condition 


SAS only 

LAMS 

SAS only 

LAMS 

SAS only 

LAMS 

248 (pilot) 

0. 0212 g 

0. 0124 g 

0. 0275 g 

0. 0166 g 

0. 0111 g 

0. 0042 g 

738 

0. 0147 

0. 0127 

0. 0217 

0. 0104 

0. 0055 

0. 003 5 

1236 (aft fuselage) 

0. 0268 

0. 0247 

0. 0403 

0. 0394 

0. 0100 

0. 0074 


The mode suppression performance described in the preceding para- 
graphs definitely shows the feasibility of demonstrating mode suppression 
techniques on the YF-12A aircraft. The canard is shown to be effective for 
mode suppression at the pilot's station, although it is relatively ineffective 
at the aft end of the fuselage. However, mode suppression at the aft fuselage 
station can be obtained through the use of inboard and outboard elevons. 

Tables XXI and XXII show the rms accelerations obtained with the canard- 
only system for various force producer input disturbances. At the super- 
sonic condition the input disturbance was introduced into the inboard elevon 
series servo to effect displacement of both the inboard and outboard elevons. 
Performance is comparable at the two flight conditions. Reduction in rms 
acceleration ranges from 73 percent at the pilot's station for supersonic 
cruise to 25 percent reduction at station 1236 for the subsonic condition. 

RMS canard surface displacements and rates were computed at the sub- 
sonic, heavyweight condition for both gust and force producer disturbance 
inputs. Table XXIII shows the computed rms canard displacements and 
rates for a canard having a total area of 3. 6 sq. ft. The canard was sized 
using aerodynamic data derived from the USAF DATCOM Handbook and is 
similar in size to the XB-70 vanes (ref. 3). 
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4.5 
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Frequency, Hz 

Figure 24. Mean Square Acceleration at FS 738 for a 1-ft/sec Turbulence 
Input With and Without LAMS System D -- F/C 1 






Figure 25. Mean Square Acceleration at FS 738 for a 1-ft/ sec Turbulence 
Input With and Without LAMS System D -- F/C 2 
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Figure 26. Mean Square Acceleration at FS 1236 for a 1-ft/sec Turbulence 
Input With and Without LAMS System D -- F/C 1 
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Figure 27. Mean Square Acceleration at FS 1236 for a 1-ft/sec Turbulence 
Input With and Without LAMS System D -- F/C 2 



TABLE XXI. - RMS ACCELERATION RESPONSE TO A 
1-RAD FORCE PRODUCER INPUT DISTURBANCE - 
SYSTEM D - F/C 1 - TRANSFORMED 18TH- 
ORDER MODEL 


Fuselage station 

Input to 
canards 

Input to inboard 
elevons 

Input to outboard 
elevons 

SAS 

only 

LAMS 

SAS 

only 

LAMS 

SAS 

only 

LAMS 

248 (pilot) 

6. 52g 

4. 45 g 

1.27g 

0. 430g 

1. 39 g 

0. 427g 

738 

1. 79 

0. 967 

0. 338 

0.227 

0. 355 

0.258 

1236 (aft fuselage) 

2. 91 

1. 53 

1. 06 

0. 791 

1. 47 

1. 07 


TABLE XXII. - RMS ACCELERATION RESPONSE TO A 
1-RAD FORCE PRODUCER INPUT DISTURBANCE - 
SYSTEM D - F/C 2 - TRANSFORMED 18TH- 
ORDER MODEL 


Fuselage 

station 

Input to canards 

Input to inboard/ 
outboard elevons 

SAS only 

LAMS 

SAS only 

LAMS 

248 

0. 62 1 g 

0. 391 g 

0. 52 6 g 

0. 143 g 

738 

0. 220 

0. Ill 

0. 222 

0. 108 

1236 

0. 300 

0. 148 

0. 376 

0. 237 


Based on the data in Table XXIII a 3a gust (12-ft/sec) will result in a 3a 
rms canard deflection of 6. 5 deg and a surface rate of 199 deg/sec. The 
total force produced by the canard in this gust environment is 985 lb. This 
force level exceeds the design goal of 600 lb. Whether or not this force level 
is acceptable can only be determined by an internal loads analysis. 

The maximum allowable forces determine the maximum magnitude of the 
test input allowed through a force producer. Inboard and outboard elevon test 
inputs were sized to produce the same rms acceleration at station 248 as a 3a 
gust input. The equivalent test input for the inboard elevon was computed to 
be 11.5 deg and for the outboard elevon. 10. 5 deg. 
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TABLE XXIII. - RMS CANARD SURFACE DISPLACEMENT 
AND RATES FOR GIVEN DISTURBANCE INPUTS - 
SYSTEM D - F/C 1 - TRANSFORMED 
18TH-ORDER MODEL 


Disturbance 

input 

Canard surface 
displacement, 
deg 

Canard surface 
rate, 
deg / sec 

1 -ft/ sec turbu- 
lence input 

0. 63 

8. 33 

1-rad input to in- 
board elevons 

32. 0 

540. 0 

1-rad input to out- 
board elevons 

27. 4 

535. 0 

1-rad input to 
canards 

1. 93 

Not 

available 


System E - Canard Plus Outboard Elevons 

Table XXIV compares roots of the canard-only system and the canard- 
plus-outboard elevon system. This comparison shows an improvement pri- 
marily in the damping of the second and fifth modes. 

Table XXV shows the stability margins obtained with the transformed 
18th-order model at the three flight conditions. At F/C 2 the stability was 
computed with the canard loop gain increased by a factor of 10 and the out- 
board elevon loop gain increased by a factor of 6. 

Table XXVI shows mode suppression performance obtained with this 
system for a turbulence input. The corresponding PSD plots are shown in 
Figures 28 through 33. These rms acceleration levels indicate mode sup- 
pression performance essentially identical with that obtained using the canard 
alone. The difference between the two systems shows up at station 1236, 
especially at the supersonic condition. The improved mode suppression can 
be seen better by comparing PSD plots (see Figures 27 and 33). The reduc- 
tion in rms acceleration at station 1236 is 43 percent versus 26 percent ob- 
tained with canard alone. It is expected that an outboard elevon system de- 
signed especially for use with a canard controller would yield further im- 
provements. 

The reduction in rms acceleration for a 1-rad white noise input applied 
to a force producer actuator is shown in Tables XXVII and XXVIII. These 
data show similar trends in performance as did the turbulence input responses. 
The only indicated improvements in rms accelerations over the canard-only 
system are at the aft fuselage station (i. e. , 1236). This is especially evident 
at the supersonic condition, where an additional 35-percent reduction in rms 
acceleration was obtained. 
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TABLE XXIV. - COMPARISON OF ROOTS OBTAINED 
WITH SYSTEM E - F/C 1 - SIMPLIFIED MODEL 


Mode 

Canard only 

Canard plus 
outboard elevon 

Damping 

ratio 

Frequency, 
rad/ sec 

Damping 

ratio 

Frequency, 
rad/ sec 

Short-period 

0. 59 

4. 24 

0. 59 

4. 24 

1st mode 

0. 202 

9. 89 

0. 196 

9. 96 

2nd mode 

0. 059 

16. 3 

0. 132 

16. 1 

3rd mode 

0. 141 

25. 9 

0. 157 

25. 7 

4th mode 

0. 059 

43. 2 

0. 077 

43. 3 

5th mode 

0. 044 

49.4 

0. 070 

49. 1 

6th mode 

0. 067 

59. 8 

0. 068 

59. 8 

7th mode 

0. 038 

73. 9 

0. 038 

73. 9 

8th mode 

0. 032 

86. 7 

0. 033 

86. 7 


TABLE XXV. - STABILITY MARGINS - SYSTEM E - 
TRANSFORMED 18TH-ORDER MODEL 


Flight 

condition 

Canard loop 
closed 

Outboard elevon loop 
closed 

Gain 

margin 

Phase 

margin, deg 

Gain 

margin, dB 

Gain 

margin, deg 

1 

— 

52 

18 

90 

1A 

— 

— 

— 

90 

2 

— 

120 

11 

72 
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Figure 28. Mean Square Acceleration at FS 248 for a 1-ft/sec Turbulence 
Input With and Without LAMS System E -- F/C 1 
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Figure 29. Mean Square Acceleration at FS 738 for a 1-ft/sec Turbulence 
Input With and Without LAMS System E -- F/C 1 
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Figure 30. Mean Square Acceleration at FS 1236 for a 1-ft/sec Turbulence 
Input With and Without LAMS System E -- F/C 1 
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Figure 31. Mean Square Acceleration at FS 248 for a 1-ft/sec Turbulence 
Input With and Without LAMS System E -- F/C 2 
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Figure 32. Mean Square Acceleration at FS 738 for a 1-ft/sec Turbulence 
Input With and Without LAMS System E -- F/C 2 
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Figure 33. Mean Square Acceleration at FS 1236 for a 1 -ft/ sec Turbulence 
Input With and Without LAMS System E -- F/C 2 



TABLE XXVI. - MODE SUPPRESSION PERFORMANCE FOR 
A 1 -FT/SEC TURBULENCE INPUT DISTURBANCE - 

SYSTEM E 


Fuselage 

station 

F/C 1, high subsonic 
heavyweight condition 3 

F/C 2, supersonic 
cruise condition 13 

SAS only 

LAMS 

SAS only 

LAMS 

248 (pilot) 

0. 0211 g 

0. 0125 g 

0. Olllg 

0. 00402g 

738 

0. 0148 

0. 0127 

0. 0055 

0. 00345 

1236 (aft fuselage) 

0. 0253 

0. 0234 

0. 0100 

0. 00570 


Simplified model 
^Transformed 18th-order model 


TABLE XXVn. - RMS ACCELERATION RESPONSE 
TO 1-RAD FORCE PRODUCER INPUT 
DISTURBANCE - SYSTEM E - F/C 1 - 
SIMPLIFIED MODEL 


Fuselage 

station 

Input to canards 


Input to inboard / 
outboard elevons 

SAS only 

LAMS 

SAS only 

LAMS 

248 

6. 52 g 

4. 44 g 

2. 156 g 

0. 7003 g 

738 

1. 79 

0. 958 

0. 576 

0. 3036 

1236 

2. 90 

1. 53 

1. S58 

0. 876 


TABLE XXVHI. - RMS ACCELERATION RESPONSE 
TO 1-RAD FORCE PRODUCER INPUT 
DISTURBANCE - SYSTEM E - F/C 2 - 
TRANSFORMED 18TH-ORDER MODEL 


Fuselage 

station 

Input to canards 

Input to inboard/ 
outboard elevons 

SAS only 

LAMS 

SAS only 

LAMS 

248 

0. 621 g 

0. 3876 g 

0. 526 g 

0. 1249 g 

738 

0. 220 

0. 1076 

0. 222 

0. 0967 

1236 

0. 300 

0. 1461 

0. 376 

0. 1574 
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Surface motion was not computed for this system configuration. It was 
assumed the surface motion computed for the canard-only system would not 
change appreciably for this system. Surface motion for the outboard elevon 
system was computed for System A, the inboard/outboard elevon system. 
These data were also considered to be representative of those for the canard- 
plus-outboard elevon system. 

Handling Qualities 

A brief analysis was performed on the analog computer to evaluate the 
effects of the LAMS, control on aircraft handling qualities. The analysis was 
performed on the analog computer using only the rigid-body pitch-and-plunge 
degrees of freedom plus the first three aeroelastic modes. The handling 
qualities were evaluated at F/C 1. 

The YF-12A primary control system uses a mechanical linkage between 
the control stick and the inboard elevon actuator. This configuration pro- 
vides a surface position command in proportion to stick position, but does not 
augment the aircraft's handling qualities. The short-period dynamics and 
static stability are augmented by the pitch SAS, however. 

A LAMS system can affect the handling qualities both in the transient 
response behavior and in the steady-state stick effectiveness. Providing 
mode suppression actually improved the transient response due to the damp- 
ing augmentation provided for the lower-frequency structural modes. This 
was true for all LAMS system configurations. The short-period frequency 
and damping were not degraded by the LAMS systems because of the design 
constraints imposed during the LAMS system synthesis. For the canard 
system, however, it was necessary to add a pilot stick input to the actuator. 
This input was necessary because the canard controller used an acceleration 
and rate feedback which acted to oppose pilot commands. This caused large 
deflections and tended to reduce the steady-state stick effectiveness. The 
addition of a pilot input to the canard actuator restored the effectiveness and 
reduced the canard displacement during a pilot command. Reduction of the 
canard displacement in turn reduced the loads imposed on the aircraft for- 
ward fuselage during pilot commands. 

Nonlinearities 

The analog computer simulation used for evaluating handling qualities 
included actuator rate and displacement limits. Generally, mode suppression 
performance is degraded by any limiting action of the actuator. The degrada- 
tion is caused by an effective reduction in force output relative to the actuator 
input command. This was not the case for the mass-reaction system. When- 
ever the mass would strike its displacement limit, its force output would be 
momentarily increased, thereby enhancing the mode suppression performance. 

Nonlinearities, such as actuator hysteresis, were not evaluated, per se, 
in the study. However, previous studies established (e. g. , see ref. 4) that 
nonlinearities such as hysteresis can result in a significant degradation in 
LAMS system performance. Hysteresis can have a significant effect because 
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surface motion required for mode suppression in the presence of a gust is 
generally small. For example, typical values of hysteresis for the inboard 
and outboard elevon actuators (and linkage) are 0. 15 deg and 0.25 deg, 
respectively. Table IX shows the rms surface displacement for System A 
to be only 0. 32 deg for the inboard elevon and 0. 67 deg for the outboard ele- 
von for a 12-ft/sec gust disturbance. Certainly, hysteresis would have a 
significant effect on the performance of System A. It is expected that by 
using a small, lightweight, dedicated surface, such as the canard vanes, the 
hysteresis effects on mode suppression performance could be minimized. 


CONCLUSIONS 


The study has shown that the application of the LAMS technology to the 
YF-12A achieves significant improvement in vehicle rms structural load 
factor for both gust inputs and for force producer inputs. Therefore, the 
YF-12A aircraft can be an effective research tool for expanding LAMS tech- 
nology and for developing improved application techniques. 

The canard vanes were found to be the single most effective force pro- 
ducer for reducing structural mode contributions to local accelerations along 
the fuselage. Use of a dedicated force producer, such as the canard vanes, 
also offers the advantage of allowing the design to be optimized for the mode 
suppression function. 

Of the five candidate systems evaluated, the one using the canard vanes 
plus the outboard elevons offers the greatest potential for continued develop- 
ment of mode suppression. There are no major technical obstacles, exist- 
ing or anticipated, that would preclude implementation of the system on the 
YF-12A. 
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APPENDIX A 
BASIC VEHICLE DATA 


The basic vehicle description used in the LAMS YF-12A study is pre- 
sented in this appendix. This description includes the symmetric structural, 
aerodynamic, and gust models for the four flight conditions that were studied. 
The choice of flight conditions was based partly on what was considered to be 
critical areas of the flight envelope from a loads and dynamics standpoint and 
also on the availability of data during the initial part of the study. The flight 
conditions studied are listed below: 


Flight 

condition 

Weight, 

lb 

Altitude, 

ft 

Mach no. 

1 

124 271 

25 000 

0. 95 

1A 

68 693 

25 000 

0. 95 

2 

90 703 

Cruise 

condition -- 

3 

68 693 

22 000 

0. 68 


NOMENCLATURE, SIGN CONVENTIONS AND 
COORDINATE SYSTEMS 


This section contains the nomenclature for the material presented in this 
appendix, as well as all sign conventions and coordinate systems used to 
describe the vehicle. 


Nomenclature 


Matrix notation. - 

[ ] rectangular or square matrix 

L 1 diagonal matrix 

{ } column matrix 

L J row matrix 

T 

[ ] matrix transpose 

Li 1 unit diagonal matrix 
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convolution operator; e. g. , 


a 



C l (t_t/) 

a 


ik. 

dt' 


dt 


/ 


J 


Coordinates and forces. - 


displacement of reference point 


rotation of reference point 


z . 
1 


e. 


displacement of i^ coordinate 
th 

rotation of i coordinate 


1 


6 r 6 2 


6 


3 


6 

m 


r. 

l 


[D z ] 


[D e ] 


Z 

O 

M 

o 

Z. 


rotation of inner and outer elevons with respect 
to hingeline, rad 

rotation of canard surface, rad 

displacement of mass of mass-reaction device 

i th generalized displacement 

transformation of deflections at structural nodes 
to aerodynamic nodes 

transformation of deflections at structural nodes 
into thetas at aerodynamic nodes 

force at reference point 

pitching moment at reference point 

force at i^ coordinate 


M. 

i 

Q, 

A LW’ a nac 


pitching moment at reference point 
i** 1 generalized force 

aerodynamic influence coefficient matrix 
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A Z X ' A Z 2 j A 0 1 ' 

An > A , A , 

0 2 g 2 

A . , A* , A. , 
6 1 6 2 6 3 


A-7 , A 


0 


generalized aerodynamic matrices 


[M], [M] 
[K], [K] 


generalized mass -reaction force matrix 

mass and generalized mass matrices 
stiffness and generalized stiffness matrices 


[D ( ] 




[T] 


{w} 


[w g 5 

Parameters . 

V 

P 

q 

b_ 


«i 


co - 
1 


transformation of relative degrees of freedom 
to absolute degrees of freedom 

matrix of mode shapes used to compute 
generalized forces 

[T] = [D.] [V ] 

1 m J 

eigenvector or mode shape 

downwash on structure, w. = V0. + z. 

i ii 

gust downwash 


free-stream velocity, in/ sec 

2 4 

density of air, lb-sec /in. 

2 

dynamic pressure = 1/2 pV 

semiroot chord, in. 

semispan of wing, in. 

damping coefficient 

i^ elastic mode frequency, rad / sec 

X. = w . 2 

l i 
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s 

s = iu = V-l u 

$ (u) 

power spectral density 

2 

a u 

2 , 2 

mean square gust velocity, ft /sec 

L 

scale of turbulence, ft 

V 

o 

velocity, ft/ sec 

u 

frequency, rad/ sec 

a., b. 
i i 

coefficients for indicial lag, i = 1, 2 


Subscript. - 

abs denotes quantity in absolute coordinate system 

rel denotes quantity in relative coordinate system 

i, j integers 

Sign Conventions and Coordinate Systems 

The basic sign conventions for forces, moments, deflections, and rota- 
tions for the symmetric analysis are shown in Figure Al. 

Two basic coordinate systems are used in presenting the vehicle dyna- 
mics. Both systems have the same axis orientation as in Figure Al. The 
first system is the "absolute" coordinate system and is defined by the motion 
of the vehicle as seen by an observer stationed on the ground. The other 
system is the "relative" coordinate system and is defined by the motion of 
the vehicle as seen by an observer stationed at some reference point on the 
vehicle, which for this study is FS 900, WL 100, and WS 32. 

These coordinate systems can be related to one another, and this rela- 
tionship will be used in the remaining sections of this appendix in defining 
the basic data. 

Transformations for displacements and rotations . - The equations 
governing the relationship between the absolute and relative displacements 
and rotations are: 

z = z 

°abs °rel 

0 =0 

°abs °rel 


(Ala) 

(Alb) 
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z, z 

Figure Al. Sign Convention 


z. = z. + z + x. 0 (Ale) 

1 abs x rel °rel 1 °rel 


0. = 0. +0 (Aid) 

l abs x rel °rel 

til 

where 5L is the streamwise distance from the i l coordinate to FS 900, with 
x. positive for coordinates aft of FS 900. 

Since the vehicle description will be given in matrix form, the above 
transformation can be written as: 



75 



Appendix A 


Transformation for forces and moments . - The equations relating the 
forces and moments in the relative and absolute coordinate systems are: 


o 


rel 


= Z„ 


abs 


n 

I 

i=l 


Z. 


abs 


M 


o , 
rel 

Z. 


x rel 

M. 


M 

°abs 


+ 


I 


i=l 




x .Z. 

1 1 abs 


Z. 

X abs 

M. 


x rel X abs 

The above equations can be written in matrix form as 


( \ 


. 




r „ \ 

z 


1 0 1 

...10... 

0 


z 

o 






o 

M 


0 1 x, 

X 1 

• • • u „ • • • 

1 


M 

o 


1 

n 



o 

Z 1 

• 


. 0 

• • 




z i 

• 

• 

< • 

> 

• • 

• • 

T 


< 

• 

• 

z 



_L 



z 

n 






n 

Mi 

■ 


• • 

• • 




M 1 

• 

• 

M 


• • 

• • 




M 

V. n J 

rel 





k n J 


(A2a) 


(A2b) 


(A2c) 

(A2d) 


(A2) 


If one examines the matrices in transformations (Al) and (A2), one notes 
that the matrix in (A2) is just the transpose of the matrix in (Al). 

The actual coordinates used in describing the vehicle dynamics are all 
displacements, except for d Q ; hence, there are no local moments or 

rotations 0^ . For this system, transformation (Al) becomes 
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r 

Z 


1 

0 

o 

e 


0 

1 

o 

z. 


1 X 

1 

• 

> = 

1 

• 


• 

• 


• • 

z 


1 X 

n 

V . 

abs 

n 


0 . 




( 




z 

o 

e 

0 
z 

1 


z 

n 





( * 1 
o 




9 

o 

\ 



z 

1 

/ = 

D. 

l 



rel 



z 

n 


(A3) 


.STRUCTURAL MODEL 


The structural model was obtained by using a two-dimensional finite 
element representation, where structural members were lumped together 
in order to keep the size of the model small but still provide a good descrip- 
tion of the dynamic behavior of the vehicle. The resulting model was then 
transformed as described below to obtain a description of the vehicle in terms 
of generalized forces and displacements. 


Structural Node Locations 

The stiffness matrix for the vehicle is defined in the relative coordinate 
system and contains 66 degrees of freedom, including the reference pitch and 
displacement degrees of freedom 9 Q and z Q . All of these degrees of free- 
dom are displacements, except for Q 0 . Their locations are listed in Table 
A1 and are also shown in Figure A2. 


Stiffness and Mass Representation 


The stiffness matrix was determined using a two-dimensional finite 
element representation, and symmetric boundary conditions were imposed. 
This matrix was transformed into the relative coordinate system by imposing 
a pitch and plunge constraint at FS 9 00, WS 32. The form of this matrix is 


r 

z 

o 


M 

o 


z i 

> = 

.V 

rel 




— 




- - 


r \ 

0 . 

0 


Z 




Z 





o 




o 

0 . 


0 


6 




e 





0 




o 

0 




Z 1 

> = 

K 

( 



Constrained 








• 

K, lb/in 



• 




• 

0 




z 




z„ 


— 




l n J 

rel 

- - 


n J 


(A4) 
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TABLE Al. - STRUCTURAL NODE LOCATIONS 


Node 

FS 

WS 

Node 

FS 

WS 

Node 

FS 

WS 

1 

900 

32 

23 

1034 

72 

45 

1212 

133 

2 

900 

32 

24 

1130 

72 

46 

795 

207 

3 

160 

32 

25 

1178 

72 

47 

914 

207 

4 

248 

32 

26 

1212 

72 

48 

1034 

213 

5 

388 

32 

27 

1236 

72 

49 

1130 

213 

6 

555 

32 

28 

1262 

72 

50 

1178 

213 

7 

640 

32 

29 

795 

100 

51 

1212 

207 

8 

738 

32 

30 

882 

100 

52 

1212 

213 

9 

795 

32 

31 

914 

100 

53 

1034 

240 

10 

954 

32 

32 

954 

100 

54 

1082 

240 

11 

1034 

32 

33 

1034 

100 

55 

1130 

240 

12 

1130 

32 

34 

1130 

100 

56 

1178 

240 

13 

1178 

32 

35 

1178 

100 

57 

1212 

240 

14 

1212 

32 

36 

1212 

100 

58 

1082 

267 

15 

1288 

0 

37 

1236 

100 

59 

1130 

267 

16 

1236 

32 

38 

1262 

100 

60 

1178 

267 

17 

1262 

32 

39 

795 

133 

61 

1212 

267 

18 

738 

72 

40 

914 

133 

62 

1130 

294 

19 

795 

72 

41 

1034 

127 

63 

1178 

294 

20 

882 

72 

42 

lieu 

127 

64 

1212 

294 

21 

914 

72 

43 

1178 

127 

65 

1178 

320 

22 

954 

72 

44 

1212 

127 

66 

1212 

320 




Appendix A 

where Z and M are the elastic force and moment at FS 900, which must be 
o o 

zero for equilibrium. 

The inertia description of the vehicle consists of a matrix of lumped 
weights at the node points. This description is in the absolute coordinate 
system and hence must be transformed. Thus, the description of the mass 
in the relative coordinate system becomes: 

[M] rei = J [D i ]T LW1 [D i ] (A5) 

where 

p 

g = acceleration due to gravity, in. /sec 

LWl = diagonal matrix of weights in lb for some specified weight 
condition 


Determination of Generalized Forces and Displacements 


Due to the lumped parameter representation of the vehicle, the general- 
ized forces can be easily determined. One first finds the nonzero eigenvalues 
and their corresponding eigenvectors for the following vibration problem: 


where 

i 


\ CM] - [K] 




(A6) 



.th 

1 


natural frequency squared of the elastic vehicle 


V\ = associated eigenvector for V 

Let [VI be the matrix described below: 
m 


“■ — 






1.0 0 
0 1.0 

0 


V 

m 

II 

• • 

• • 

0 0 


- 


- 



' 1 


\ V 

> 

n 





( A 7 ) 


From the above matrix, the generalized forces and displacements can be 
defined for the vibration problem as: 
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n+2 

z. = y~ V r. 
1 Z_ m. . ] 

3=1 13 

n+2 

Q. = Y V F. 

i m. . i 

• i ii J 

3=1 J 


where Q. are the generalized forces and r^ are the generalized displacements. 

The above transformation can be carried out by use of [V ] and its trans- 

J m 

pose. The form of the generalized mass and stiffness matrices becomes 


[K] = [v m ] T Dec] [v m ] (as) 

[M] = [V m ] T Dvt] [V m ] (A9) 


The same procedure is used for converting the aerodynamic matrices, 
since the airloads are lumped at many stations and, hence, is a good approxi- 
mation to the integral form of the generalized forces. 

Sixteen modes, corresponding to the 16 lowest-frequency symmetric 
modes of the free -free vehicle, were used to modalize the data for this study. 
This representation allows one to work in the 0- to 25-Hz range for dynamic 
calculations. 


Representation of Structural Damping 


The structural damping of the vehicle is assumed to be proportional to 
the displacement and in phase with the velocity. This representation is 
usually carried out by multiplying [K] by a diagonal matrix of the form 
[l + ig l, where g. is the damping coefficient for the j** 1 mode. In this study, 
3 3 


it was advantageous to keep all matrices real; hence, an equivalent viscous 
damping formulation was used. Since the natural frequencies of the system 
with and without aerodynamics are more or less constant, the equivalent 
viscous form of the damping can be approximated as 



x. K.. 

3 3 33 


CJ. 

3 


d. . = 0 for i ^ i 
13 J 
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where 

k.. 

33 


co. 

3 


The damped vibration equation is then written as: 

+ [K] {r} = {0} (A10) 


0 

0 


0 


{r} = generalized coordinate or displacement vector. 


where 


[M] £r} + [D] {r} 


0 0 . . . 


[D] 


d.. 

IX 


0 0 L 


f Vi 

= j generalized stiffness 
= j** 1 vibration natural frequency, rad/ sec 

= damping coefficient assumed to be 0. 02 for the first two modes 
and 0.05 for the rest 


AERODYNAMIC MODEL 

The aerodynamics used in this study were formulated in the time domain 
instead of the frequency domain. This formulation was needed for the control 
theory problem discussed in Appendix B. 


Subsonic Aerodynamic Model 

For the subsonic flight conditions, the aerodynamics for the wing and 
forebody were based on the subsonic kernel function or lifting surface theory 
(ref. 3), and used in conjunction with a uniform indicial lag. The airloads on 
the nacelle and those on the vehicle due to elevon rotation were determined 
from wind-tunnel measurements and were used without a lag function. 

A steady-state aerodynamic influence coefficient matrix for the wing and 
forebody was computed, and adjustments to the airload distributions, based 
on pressure measurement tests, were made. This matrix is used to repre- 
sent the airloads at the structural degrees of freedom due to the displace- 
ments and slopes of the structure at the nodes. 

To account for the lag in the buildup oflift due to the motion of the struc- 

-b.t 

ture, a uniform lag of the form 1 - a^ e was used. The values of a^ and 

b^ were estimated from ref. 7. The effect of this lag was checked when a 
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flutter analysis of the vehicle was performed using both the uniform lag aero- 
dynamics and the unsteady kernel function aerodynamics. That analysis 
showed that the flutter results for the vehicle were basically the same for 
both representations. 


The form of the indicial aerodynamics is outlined below: 



where 


(All) 


w(t) 

V 

a lw ( ‘> 


represents ( the convolution operator 
semiroot chord lengths traveled 

nondime nsional downwash at the structural degrees of freedom 



The above equation can be transformed into the frequency domain and has 
the following form: 




{z i 3 


+ 


abs 


4tt£ 2 q 


^ a nac 3 


Ki 



rel 


(A 12) 


where 

s = iu 

j = represents a transformation of the deflections, z^, at the 

L structural modes into slopes at the collocation points used 

in the formation of the aerodynamic influence coefficient 
matrix [Q(iu)], i = 1 for wing forebody, i = 2 for nacelle. 
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similar to 



except it relates the deflections at the struc- 


l L J 

tural modes to the deflections at the collocation points. 


Since the above formulation is in the absolute coordinate system, the 

transformation [d^] can be applied to transform into the relative system. 

This transformation and the reduction to generalized forces and displacements 
will be performed simultaneously by using the following transformation: 

[ T ] ■ [ D i] [ V m] 

Under the above transformation, one has 


!A(s) V = 4n^ 2 q [t 1 
aero I re ^ L J 


+ '“NAC 

or for fixed q and V 


1 - 


a-^s 


s+b. 


IV_ 

b 


D, 


+ -S- 
V 


o i 


Dr 


[ a lw]( 


D, 


V 


Dr 


(A1 3) 


[t] {r} 


rel 


}A(s) \ 
aeroj re ^ 


1 - 


ii 


■*> 1 1^1 


<\ +s > 


A Z, I ) + | A 0, 


+ s 


{r} 


rel 
(A 14) 


The aerodynamics due to elevon rotation were needed in the study. The 
forces due to the inboard and outboard elevons were calculated separately and 
were not lagged. These forces were determined from wind tunnel test data 
and include the induced effects on the wing itself, as well as surface forces. 
The form of these forces is given below: 



(A15) 


(A16) 
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where 

6 ^ = surface angle of inboard elevon relative to wing 

62 = surface angle of outboard elevon relative to wing 

df, $2 = positive for elevon trailing edge down 

By use of the transformation, T, one obtains the generalized forces due 
to elevon rotation, for fixed q, as 


K! 

; = q[T]T i L sJ 

; 6 i = : 

Kh 

(A17) 


rel ' 



h; 

L ■ * mT N 

M 

\h 

(A 18) 


The gust aerodynamics for the study presented a problem since they had 
to be formulated as indicial-type aerodynamics. An interim uniform lag 
function of the form 


1 


e 



t 


was used, and the value of b 2 was based on results presented in refs. 5, 7, 
and 8 . 

The interim gust aerodynamics was defined to be: 



,w 

( 


g 


V 

+ [a nac ] ” 

V 


f dt 

A 

dt 



(A 19) 


or using the transformation, T, one obtains the forces in relative system for 
fixed q and v, as 

T 

S 


fQ gust (s) Vel 


1 - 


s+ b r 


V\ 

,TT* 

o 


w_ + )A 


Si l S 


w 


g 2 8 


(A20) 
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Supersonic Aerodynamic Model 


For the supersonic flight condition, i. e. , the cruise condition, the aero- 
dynamics were formulated using piston theory due to the high Mach number. 
Because of the Mach number considered, no aerodynamic lag was used on 
the aerodynamics due to the motion of the structure. A lag on the gust aero- 
dynamics was used in the gust analysis of this vehicle and was retained for 
this study. 


The form of the aerodynamics is basically the same as that used for the 
subsonic conditions, except that the lag term is deleted in all but the gust 
formation. Thus, one has: 


Khbs = 
(Slabs = 

{ Z B ust labs = 



(A21) 

(A22) 

(A23) 

(A 24) 


By applying the transformation, T, the above forces are transformed to 
the relative system and become generalized forces and displacement for 
fixed q, V 


Kft 

(s)\ = - 

J rel 

' { + V } k^rel 

(A25) 

K 

1 , ■ 

{\} 6 i 


(A26) 

K 

1 ■ { 
rel 

: A 6 2 } 5 2 


(A27) 
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(A28) 
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EQUATIONS GOVERNING THE ELASTIC VEHICLE 


This section contains the equations governing the elastic model, and 
describes the internal loads calculation and gust models. 


Equations Governing the Vehicle Motion 

The equations of motion for the elastic vehicle without the baseline pitch 
SAS are given as: 


's [M] + 


1 - 


a^s 


s + 


■ 

/ 



r 1 

\ 


\ 

1 V\ 

r 

A z 

+ 

A e 

r s 

A z t) 

+ CD] 

' °“l 

0 1 

\ 

1 j 


L 1 

1 

2 

/ 


+ 


+ [K] 


{r} = (F(s)} 


(A29) 


where s = iu. 


The baseline SAS can be included in the above equation by incorporating 
the transfer function defined from the SAS block diagram (see Figure 2 of 
the main text). 

The different types of forcing function that could appear on the right-hand 
side of the above equation are given below in the frequency domain form. 


Elevon generalized force description. - 


{F(s)} = Ja 5 J 6. (s) i = 1, 


(A30) 


where S.(s) is the rotational coordinates for the elevons and may be a function 
of s. 


Gust generalized force description . - 


{F(s)} 


1 - 


o 


A g 1 rg (s)+ i A ^2^ w s (!i, 


(A31) 


Canard generalized force description. - 


[F(s)5 = qS C It] « 3 = (A ) 63 

6 3 


(A32 ) 


87 



Appendix A 


where 

t. 

i 


th T 

= element of the i row of the fourth column of [T] , 

th 

i. e. , displacement of the 1 mode at the location of 
the canard, FS 248 

= exposed area of one side of canard 


C T = lift slope coefficient for canard, rad 
«3 

Mass reaction device generalized force description. - 


where 


{F( s ) } = s 2 M (t- } 6 
J L 1 m 



M = mass of the mass-reaction device 
6^ = displacement of M 

t. = identical to L in the canard description 


(A33) 


It should be noted that various combinations of the above forcing functions 
can be used, with the exception of the case involving both the canard and mass- 
reaction devices. These devices are considered mutually exclusive of one 
another. 


Description of the Gust Model 


The preliminary LAMS longitudinal gust model contains the w ind filter, 
associated lift growth functions, and gust force coefficients. There are no 
penetration effects included in the present model. The input to the wind 
filter, r|, is a unity rms "white noise" signal. The gust model is shown in 
Figure A3. 


The wind filter is derived from the Press-Meadows power spectral den- 
sity random gust model (ref. 7), as follows: 


where 


o 


u 


2 


$ (u) 



1 Lj 

1+3(L/V ) 2 u 2 
o 

V 

\ ol 

_[l + (L/V o ) 2 u 2 ] 2 


ft 2 /sec 2 
rad/ sec 


2 2 

mean-square gust velocity, ft / sec 


(A34) 


L = scale of turbulence, ft 
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Wind filter 


Note: 


{F } and 


{F 3 contain 
g 2 



Figure A3. Preliminary Longitudinal Gust Model 


V q = velocity, ft/ sec 
w = frequency, rad/ sec 


Using the fact that §(to), the Press-Meadows power spectral density 

2 

model, is equal to w (iio) , a Laplace transfer function equivalent to Press- 

o 

Meadows model can be developed. Note that 


$(io) = |w (iu)| 2 = w (iu) w * (iw) (A35) 

& s s 

where the asterisk denotes the complex conjugate. Factoring $ (to) into appro- 
priate terms and replacing ito by s, the Laplace operator, it is found that 


w (s) = a 

o 


JK 

u V L 


V 


s + 


V3l7 


V 


s + 


or, substituting V q = 978. 33 ft/ sec and L = 1000 ft. 


w (s) 

o 


1. 713 a 
u 


s + 0. 5648 
(s+ 0. 9783)' 


(A36) 


(A37) 
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Transformation of the Vehicle Equations of Motion 

A change in the choice of generalized coordinates was made to eliminate 
zero eigenvalues from the unforced equations of motion, eq. (A29) with 
{F(s)3 = 0. This change in variables was made by choosing w Q = z Q + V0 q 

and 6 0 as the first two generalized coordinates instead of z Q and 0 o as 

are used in eq. (A29). 

Upon inspection of the components of eq. (A29), one notes that all of the 
aerodynamic forces and moments associated with the rigid-body motion 
(z Q , 0 O ) are in fact only functions of the rigid -body downwash, w Q . Further 

examination of the system reveals that the aerodynamic terms are the only 
terms that involve z Q , d Q and their time derivative except for the rigid - 

body equations which have terms involving z Q and 6 . These acceleration 

terms can be easily represented in the new system by noting that 



Carrying out the above change in generalized coordinates and dividing 
each equation by the appropriate generalized mass, one can rewrite eq. (A29) 
as 
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{F(s)} 



where F(s) is defined in eqs. (A30) to (A33) and the symbol " denotes 
the transformed matrices as outlined above. 
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DESCRIPTION OF SYNTHESIS TECHNIQUES 


This appendix describes two synthesis techniques applicable to the design 
of LAMS systems. These two techniques are referred to as the frequency 
response technique and the quadratic optimization technique. Both techniques 
have been thoroughly described in previous LAMS reports (see refs. 4 and 6). 
These descriptions are repeated in this appendix for the convenience of the 
reader. The description of quadratic optimization theory also includes the 
details of the theory peculiar to the YF-12A application. 

An attempt was made in the LAMS YF-12A study to integrate the two 
techniques into a single unified procedure. It was the intent in performing 
this integration to capitalize on the strong point of each approach. The 
attempted integration of these techniques is discussed in this appendix following 
the discussion of the two individual techniques. 


SYMBOLS 

z normal acceleration at an aircraft station where the acceleration is 

c 2 
to be reduced, in. / sec 

•• 2 
z aircraft motion measured by LAMS system, in. /sec 

s 

w gust velocity, ft/ sec 

S 

5 control surface deflection, rad 

K c G(s) gains and filtering of LAMS system 

T desired attenuation of normal acceleration with LAMS system 

[ ] rectangular or square matrix 

l 1 diagonal matrix 

( } column matrix 

L J row matrix 

[ 3^ matrix transpose 


92 



Appendix B 


convolution operator e. g. 


[ a(t)]* {b ( t ) } = < [ a(t-t') b(t') dt' 


generalized coordinate vector 


state vector 


control vector 
white noise 
response vector 
state-transition matrix 
control- input matrix 


disturbance- input matrix 


state-response matrix 
control-response matrix 


A A A AAA A A A A a 

M - A V A - B - V V V V A c- A m 


(See Appendix A) 


A A A 

C = A + A 


E = A„ + B 
d l 

AAA A 

T = b, A - A„ 

1 zj 


p = Lp^ . . . , P i g -f state variables associated with lift lag 
Pjg state variable associated with gust lag 

Dn n state variable associated with wind filter 


P2i» P 2 2> P23' P24’ ^25 state variables associated with the baseline SAS 
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w = z + V fl , in. / sec 
o o o o 


o 
z o 
r 3 ’ 

«1 

69 


' m 


u 


m 


pitch at vehicle reference point, rad 
plunge at vehicle reference point, in. 

, , r 1Q generalized coordinates, in. 

inboard elevon deflection, rad 

outboard elevon deflection, rad 

canard deflection, rad 

generalized coordinate associated with mass reaction device, 
control input to mass reaction device 
turbulence scale, ft 


K(t) = a x e 


-b^t 


\f/( t) = b 2 e 


•V 


V 


a^, b^, b 2 coefficients of indicial lags b = b 

u o 


V , V ' velocity of vehicle 
0 0 


6 


*; 

* i. 


Laplace operator 
scheduled gain 
th 

slope of the j mode shape at the rate gyro location 

th 

deflection at station j of the 1 mode shape 


rad 
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FREQUENCY RESPONSE TECHNIQUE 


This synthesis technique (ref. 6) is described by first presenting the 
general LAMS problem. The problem is presented in terms of the design 
objectives and design approach considerations which form a basis for the 
development of a LAMS system. Using this information, the formulas and 
conditions used in the synthesis procedure are established. 


Design Objectives and Design Approach 

The following design objectives should be used to form the basis for 
developing an effective LAMS system: 

• Reduce response amplitudes of those modes (generally rigid- 
body and lower- frequency structural modes) that are significant 
in terms of their contributions to rms stress and acceleration 
levels at aircraft stations of interest 

• Satisfy stability requirements 

• Maintain adequate rigid-body dynamics in terms of good handling 
qualities 

• Provide maximum tolerance to unpredicted variations in vehicle- 
system dynamics 

• Minimize unfavorable coupling with structural modes above 
frequency range of significant structural modes 

• Establish analytical design which permits economical and 
reliable hardware implementation 

A necessary condition for load alleviation and mode suppression is the 
ability to sense the deflections, or a derivative of the deflections, of the 
significant modes and to apply control forces to these modes. The deflection 
rate and/or acceleration of the modes can be measured by conventional sen- 
sors (such as accelerometers or gyros) suitably located on the vehicle. A 
corrective control force can be applied through the deflection of one or more 
of the aerodynamic control surfaces. Hence, a design approach was assumed 
which would concentrate on the development of a system using conventional 
sensors to measure the modal deflection rates and/or accelerations and to 
command corrective action through the aerodynamic control surfaces. 

Development of synthesis technique. - In general, the synthesis tech- 
nique used to design a LAMS system should have the following properties: 

• Facilitate design studies using a high-order mathematical 
model of the vehicle 
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• Enable use of this model in a frequency response form, thereby 
facilitating incorporation of unsteady aerodynamic effects in 
the frequency domain 

• Capability to compute required controller characteristics 
directly from specification of meaningful design criteria 

• Permit rapid and accurate evaluation of non- ideal controller 
characteristics 

The synthesis technique described on the following pages is an attempt to 
provide a design procedure with these properties. 

The technique consists of (1) specifying the control system filter, and 
(2) specifying the sensor configuration. The basic formula used in the 
synthesis technique specifies the control system filter (i. e. , compensation) 
requirements in terms of a given sensor configuration and an allocated 
attenuation requirement on the rms stress or acceleration level at a given 
aircraft station (e. g. , the pilot's station). For the sake of this discussion, 
only reductions in acceleration levels will be discussed. Clearly, the same 
techniques can be used for reducing stress levels. Additional requirements 
are specified in terms of stability constraints and in terms of restrictions 
resulting from the use of practical (in terms of mechanization) linear filters. 
These additional requirements enable the designer to specify a sensor con- 
figuration. Asa consequence, once a reduction in rms acceleration level 
has been established at the given stations, it is possible to specify a sensor- 
filter combination. The synthesis procedure requires the vehicle dynamics 
to be described in a frequency response format in which output variables of 
load accelerations are related to input variables of wind gust and control 
surface deflection. Specification of data in this form facilitates incorporation 
of unsteady aerodynamics in the frequency domain rather than in the time 
domain, a factor often not included in other techniques. 

The approach to be taken toward performance specification for purposes 
of controller synthesis is to specify the desired reduction in acceleration at 
the resonant peaks of each of the significant modes. The frequencies of these 
peaks will be determined from a frequency response plot between accelerations 
and the wind gust input, thereby including all known coupling effects. The 
required percentage reduction in amplitude at each frequency will be estimated 
from power spectral density plots of acceleration due to turbulence at selected 
vehicle stations. Knowing the rms acceleration and identifying its source in 
terms of individual modes, a design goal for reduction of resonant amplitude 
can be allocated to each mode. An important part of this allocation is defini- 
tion of modes which do not cause significant acceleration contribution, so that 
the required control action at those particular frequencies can be relaxed. 

The objective of this process is to produce an appropriate allocation of control 
effort. 

With this information, the controller can be specified (as will be shown) 
and then evaluated at other aircraft stations, and stability margins can be 
established. If all constraints are not met, the allocation of required reduc- 
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tions in amplitudes at the resonant modes can be modified and the process 
repeated until all the constraints are satisfied. This results in an iterative 
process but with an orderly and systematic means of attack. Relating peak 
accelerations to rms accelerations and specifying reductions required in terms 
of acceleration at a particular aircraft station are useful methods for expe- 
diting the control synthesis process. The final measures of control effective- 
ness are the actual rms accelerations at the various aircraft stations. 


To establish the basic formulas, consider a general block diagram of the 
control system /vehicle dynamics as shown in Figure Bl. This diagram 
illustrates the relationship (in transfer function form) between the accelera- 
tion level measured at a sensitive aircraft station, z , and the wind input. 


w . Clearly, the diagram could be expanded to include all the aircraft 

O •« 

stations, but this is not necessary for the following discussion. The z term 

s 


is the quantity to be measured by a sensor located at some vehicle station. 
The K c G(s) term is the filter to be designed for the LAMS system. Thus, it 


is the task of the synthesis process to determine what quantity will be sensed 

(z ) and to define a filter [K G(s)] that will result in a specified reduction 
s c 


in acceleration (z c ) in response to the turbulence input, w^. 


The aircraft 


model must be known in a frequency response form relating output variables 


of acceleration to input variables of wind turbulence. 



and control 


surface deflections, 5. Actuation system dynamics can either be included in 
the aircraft model or be left to be specified as part of the filter [K c G(s)]. 



Figure Bl. Vehicle /System Model General Block Diagram 
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With the design goals in terms of a reduction in acceleration at a selected 
station with augmentation to that without augmentation, one can derive an 
expression for the required filter characteristics [K c G(s)]. Let T equal the 

ratio of acceleration at a particular station with augmentation to that without 
augmentation. Acceleration occurring at this station (z ) without augmenta- 
tion is 



/ " \ 

/ .. \ 1 

' z c 

z. no 
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r g / 
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Acceleration at the sensitive station with augmentation is 
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Thus, by definition of p we obtain: 
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(Bl) 


(B2 ) 


(B3) 


This definition of p implies that P can be a complex number. The 
designer is free to pick any r, real or complex, within certain limits. 
Clearly, if the control is to reduce the amplitude of the acceleration at the 
discrete frequency points, then the absolute magnitude of r chosen must 
be less than one (i. e. , 0 £ T < 1). Selecting real values for p will result 

in lower gain values for the filter [K c G(s)l than will the selection of complex 

values. Low gain values for the filter are desirable from the standpoint of 
minimizing coupling with high-frequency structural modes. 


We can solve eq. (B3) for K G(s ) to yield: 


K c G( s ) 


1 

r i 

*s 

i 1 1 

/z c /g | 

fz s /w g \ 

6 

1 " i-r | 

l 2 0 /w gj 

iL/a 


(B4) 


This equation defines the filter characteristics [K c G(s)] in a frequency 

response form. It implies that if one picks the quantity to be measured (z ), 

s 

and a value of Y at discrete frequency points (e. g. , p^, r£, Tg corres- 
ponding to the first three structural mode resonant frequencies), then a gain 
and phase requirement can be computed for K c G(s ) at these discrete fre- 
quencies. If a filter can be designed to exhibit these gain and phase 
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characteristics, then the acceleration (z ) will be attenuated by the values 

of r chosen at these frequency points. It is sufficient to specify a filter 
requirement at only the resonant frequencies. Thus, the resulting filter will 
be simpler than the inverse of the vehicle transfer function which would 
probably be obtained if the accelerations were to be reduced by only a cancel- 
lation technique. 

The remaining task is to establish a procedure for defining the quantity 

to be measured (z ). The objective is to specify the relative modal compo- 

s 

• • 

nents of z and not the type of sensor nor the combination of sensors required 
s 

• • 

to actually provide z . The sensor combinations can be established by 

s 

• • 

examining the properties of z and comparing them with the mode shapes. 

s 

The term z is considered to be an acceleration quantity having the following 
form: 

n 

+ Z K i' ; i+2 (B5) 

i=l 

where 

z q = rigid-body plunging mode 

B q = rigid-body rotation mode 
• • 

r . = structural modes 
1 

• • 

In the discussion that follows, z is referred to as the sensor comple- 

s r 

ment. Thus, to define z corresponds to defining the relative values of the 

s 

modal deflections (K^). To do this, we examine the following constraints: 

1. A stability constraint is imposed on the system by the stability limits 
of the feedback loop shown below (see Figure Bl): 



Figure Bl. Vehicle/ System Model General Block Diagram 
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This loop has the closed-loop transfer function given by: 

K G(s ) 

H(S) = 77^ (B6) 

1 - (z / 6) K G(s ) 

b L 

The term H(s) is unstable whenever the roots of its demoninator lie in 
the right-half plane of a root locus diagram. The denominator will be unstable 
at the frequency at which 


0 -. + 0 
w s5 


0 deg 


(B7) 


and 



where 


0 _ = phase angle 
s o 


of transfer function (z / 6) 

b 



= phase angle of filter [K c G<3)] 


(B8) 


Since one cannot reduce the amplitude of the structural mode accelera- 
tions by gain stabilization, it is necessary to provide phase stabilization by 
proper selection of 0 g g + 0g. Because of uncertainties in the mathematical 

model, it is necessary to require that over the range of frequencies where 
gain stabilization is not attainable nor desirable, a phase margin constraint 
is applied. For example, if it is desired to have minimum phase margins of 
90 deg, the constraint becomes: 

90 deg £ 0 _ + 0 £ 270 deg (B9) 

' s 0 g 


At frequencies above the modes of interest (i. e. , above frequencies 
where augmentation is to be applied) it is desirable to provide gain stabiliza- 
tion to eliminate this constraint on the phase angle. 

With the above constraint on the phase angle of the sensor/filter combina- 
tion, the problem is to find a solution which, over the frequency range of 
interest, will simultaneously satisfy this constraint and the phase constraint 
imposed by the filter equation [eq. (B4)]. These two constraints are repeated 
here for reference: 


K GM - 
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(BIO) 
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90 deg <: 0 g5 + 0 g ^ 270 deg 


(BIO) 


To impose the constraints defined by eq. (BIO), we observe that 


0 + 0 g = phase angle of 


Defining 


r i 

1 ' r=r 



lV w g / 

\V 6 / 


N = 


V 



(Bll) 


0 1 = [phase of (z 1 6) - phase of (z / w )] 

X C , c g 

• • •• 

0 = [phase of (z /w ) - phase of (z / 5 )] 

Z s 6 s 

and restricting 0 g g + 0 g according to eq. (B9), we obtain a minimum attain- 
able value on F, the attenuation factor: 

r 2 1 - VN cos (0 X + 0 2 ) 


(B12 ) 


At any flight condition, N and 0 . are defined for a given vehicle station. 

Thus, if we pick a quantity to be measured (z ), we can determine the lower 

s 

bound on r (defined as F ■ ) in order to assure that the stability constraint, 
1 m in u > 

eq. (BIO), is satisfied. That is: 

F * F min = 1 - VN cos (0 X + 0 2 ) (B13) 

If we assign a value to r , say F^, at a given frequency (u^), then we 

must find a quantity to be measured (z ), such that 

s 


r i* 


r . (u>.) 

min l 


(B14) 
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or else we cannot simultaneously satisfy the stability constraint and provide 
the attenuation, r = I\. The above condition is a necessary condition for 

specifying "z , but it is not a sufficient condition, 
s 

2. To further specify the requirements for the quantity to be measured 

(z ), it is necessary to examine other design objectives and analyze the 
s 

practical constraints imposed on the filter [K c G(s)]. The following design 

objective imposes a severe constraint on the type of filters and on z : 

s 


Minimize unfavorable coupling with structural modes which lie 
above the frequency range of the modes to be controlled. To 
assure this objective is met, we impose the condition 



K G(s) 
c 


< 0. 5 


at modal frequencies above the highest significant mode. 

3. It is assumed the filter [K c G(s)] will be a linear filter made up of 

first- or second- order elements and that the poles of the filter will be stable. 
Thus, the design objective of high-frequency gain stabilization, for practical 
purposes, imposes an additional constraint on K c G(s). To meet this objec- 
tive, the gain of the filter will be decreasing with increasing frequency at a 
rate of at least 6 dB per octave over the frequency range of the significant 
modes (for this argument assume the significant modes are the first three 
modes). From this, it becomes evident that if the same amount of attenuation 
(or more) is required at the second or third mode than as is required at the 
first mode frequency, then we impose 



for i 


2, 3 


(where BG is the i mode deflection) to compensate for the loss in gain due 
to the filter [K c G(s)]. For gain stabilization of the fourth and higher modes 
(assuming three significant modes), it is generally desirable to minimize BG 
for i = 4, 5, ... 


4. A further consideration is the amount of phase or gain change with 
frequency required of the filter [K G(s)]. Sensor configurations should, in 

general, be chosen to minimize large phase or gain change requirements with 
frequency for the filter to facilitate the hardware implementation. 

5. In some applications it is an objective to minimize the coupling with 
the rigid-body mode. This can be done by placing the sensor in a position 
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where the accelerations due to flexibility are the largest with respect to those 
due to rigid-body motions. The coupling can be further reduced by minimizing 
the gain of the filter at the rigid-body frequencies. 

All of the above considerations enable one to essentially define the desired 

sensor complement, z . To summarize, the sensor complement should be 

s 

chosen to satisfy the following constraints: 


1 ) 


r . 

min 


^ r. 

i 


2 ) 



2, 3 (assuming 3 significant modes) 


3) Minimize variation in gain and phase changes with frequency 
required of the filter [K c G(s)] 

4) Minimize gain of sensor-filter combination at frequencies beyond 
significant structural modes 

Once the quantity to be measured (z ) has been defined (according to 

s 

the above rules), realistic sensor locations and combinations of sensors can 

be examined to find one which will best provide z . With the actual sensor 

s 

configuration defined, the filter requirements can be computed for that sensor 
configuration through the use of eq. (B4). Using these filter requirements, a 
practical filter can be designed. Finally, with the filter designed, the prob- 
lem can be worked in reverse to check the "goodness" of the entire system. 
Thus, compromises in filter design can be analyzed to determine their effect 
on performance. It must be kept in mind that eq. (B4) does not have to be 
satisfied at all frequencies nor is its solution the only solution to the problem. 
Once a solution has been obtained, variations can be made on it to yield a 
possible better solution. 


QUADRATIC OPTIMIZATION TECHNIQUE 

Quadratic optimization techniques are derived from the concepts of 
modern control theory. Quadratic optimization differs from the more con- 
ventional servomechanism analysis of feedback control in that it is a proce- 
dure which is applied in the time domain rather than in the frequency domain. 
Optimal control synthesis is based on the notion of minimizing an analytical 
expression used as a measure of system performance. In quadratic optimiza- 
tion this expression, called the cost function, is a linear summation of quad- 
ratic elements. Each element represents a performance parameter to be 
minimized such as pilot's acceleration. Minimizing the cost function yields 
a set of gains on a linear combination of feedbacks which makes up the con- 
troller. The feedbacks, in the ideal case, include all the state variables 
(e. g. , pitch rate, normal acceleration, etc. ). 
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To use optimal theory, it is necessary to rewrite the equations of motion 
as a set of first-order linear differential equations. This set of first-order 
equations is referred to as the standard form. Reducing the equations to 
standard form generally varies from application to application. Details of 
this reduction for the YF-12A are described in subsequent paragraphs follow- 
ing the general description of quadratic optimization. The LAMS problem is 
treated as being statistical in nature. Hence, the equations of motion and the 
response equations are recast in terms of covariance equations. The response 
equations are simply linear combinations of state variables such as the output 
of an accelerometer located at some station on the vehicle. The standard 
form for the equations of motion and the response equations is as follows: 

x(t) = F x(t) + Gj u(t) + G 2 T|(t) (B15) 

r(t) = H x(t) + D u(t) (B16) 

where x(t) is the state vector, u(t) is the control vector, r(t) is the 
response vector, and r|(t) is the white noise disturbance input. The terms 
F, Gj, G 2 > H, and D are matrices defined by the aircraft dynamics. 

To treat the problem statistically, the expected mean square values 
(covariances) for x(t) and r(t) are derived from eqs. (B15) and (B16). The 
following notation will be used: 

X(t) = E {x(t), x(t)^} = expected value of x(t) 

R(t) = E {r(t), r(t)^} = expected value of r(t) 

Then, the state vector expected value for the system with no control [u(t) = 0] 
satisfies the expression 

£(t) = F X(t) +X(t)F T + G 2 N G 2 T (B17) 


where 

E Cn(t), ri(t) T } = N 6(t - t) 

(6 is the Dirac delta function); and the response vector mean squared value 
satisfies the expression 

R(t) = H T X(t)H (B18) 

The optimal controller for a linear system is a linear combination of states, 

u(t) = K x(t) (B19) 
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The optimal control problem is that of selecting the set of gains, K, such 
that performance measure J* (the cost function) 

J* = E { r(t) T Q r(t)} (B20) 

is minimized where 

r(t) = (H + DK) x(t) (B21) 

for the controlled system. The solution is subject to the constraint 

X = (F + Gj K) X + X(F + G x K) T + G 2 N = 0 (B22) 

The weighting matrix Q in eq. (B20) is diagonal and assigns the control 

emphasis to the desired responses. 

The solution to the problem expressed by eqs. (B20) and (B22) is given 
as the solution to the two equations: 

t 

K = -(D T Q D) 1 (D T QH+G 1 T S) (B23) 

S = (F + Gj K) T S + S(F + Gj K) + (H + D K) T Q(H + D K) (B24) 

where S is a matrix of Lagrange multipliers. 

Taken together, these equations represent the matrix Riccati equation. 
For the LAMS program, eqs. (B23) and (B24) were solved by converting them 
to difference equations and then computing S and K iteratively until the gains 
converged to a steady- state value. 


Reduction to Standard Form 

The vehicle equations of motion described in Appendix A must be reduced 
to standard form in order to apply quadratic optimization. This procedure as 
it applies to the YF-12A is described in the following paragraphs. 

The equations of motion without the baseline SAS, as given in Appendix A, 
eq. A29, are given below in the frequency domain as: 
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<!M s 2 + 


1 - 


a^s 

A 

s + b. 


A Zj S + A e 


+ As + B, 


»r(s)> 


' a 6j[ 6 i + 6 2 + 


f 


s + b r 


! M w g + A g 2 w g 


+ 


A vg 0 for canard controller 

I CIO 


or 

_2 


W 


for mass reaction controller 
m 


V 


The above equation can be rewritten in the time domain as: 

M r(t) + C r ( t) + E r (t) - A K(t)* r (t) - A fl K(t)* r (t) 

Z 1 y l 


= A 


«! 5 l (t)+A 6 2 6 2 U, + A gl ^ (t) * *g + A g 2 w g 


I * 

A g n (t) for the canard controller 
c o 




or 

A • • 

A g (t) for mass reaction controller 
m m 


where 

* * 

C = A + 




A A A 

E = A _ + B 
0 1 


K(t) = a 1 e 


-bjt 


i//(t) = e 


-b 2 t 


* denotes the convolution operator 
rT = Lw o' »o' r 3 r is j 


(B25) 
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By the introduction of new state variables, eq. (B25) can be put into the 
required standard form. The state variable forms of the convolutions appear- 
ing in eq. (B25), as well as the transformation of the gust, elevon, canard, 
and mass reaction device transfer functions to the desired form, are outlined 
below. 

T 

State variable form of the convolutions . - Let p = LPj, . . . , P^gJ and P ig 

be new state variables defined as follows: 

. . - f -Mt-i) 

p = K(t)* r(t) = K(t) * r(t) = - a 1 b 1 e 1 r(T)dT (B26) 

o 


One can readily see that p is a solution of the first-order differential equation 

p = bjP - a^r (B27) 

Since k(t) *r(t) = -b^p, the terms of eq. (B25) involving the convolutions 
K(t) *r(t) and K(t) * r(t) can be written as 


-A K(t)*r(t) - A fl K(t)*r(t) = |l - A fl | = T 

z i 9 i \ z i e ij p 


(B28) 


In a similar manner, the state equation governing i^(t)*w (t) is found to be 

o 


where 


Thus 


i>19 = £ 2 p 19 ' S 2 w g (t) 
p i9 = m*" g w 


A g/ (t) * V g (H ’ A Sl P 19 


(B29) 


(B30) 


Eqs. (B27), (B28), (B29), and (B3 0) define the transformation of the con- 
volutions in eq. (B25) into the desired form. 


State variable form of the gust function . - The gust transfer function given 
in Appendix A, eq. (A36), has the form: 


w g (s> = 


/ 3V 

Vv 1 



V 

o 

VTi 

lJ 

v 1 

o 

2 


L 


T|(s) 


(B31) 
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One can easily verify that w (s) is a solution of the following second -order 

S 

differential equation: 


w (t) + 2 

o 


V 

T 


HvH'elVV? 


V 


r|(t) + 


V3 L 


T|(t) 


(B32) 


By introducing a new state variable, PgQ, eq. (B32) can be placed in the 
desired standard form as: 

V /3V - 

w g (t) = " 2 “l w g (t) + p 2o + V r 2 ^ (t) 


(B33) 


where p 2 q satisfies 


20 


V 

_o 

L 


V 


w g (t) + — q(t) 





(B34) 


State variable form of the elevon equations of motion . - The transfer 
function for the inboard elevon is shown in Figure B2. Upon inspection, one 
finds the governing differential equation for 6 ^(t ) to be: 

6^ = -50. 5 6 1 - 1568 6 X + 1568 Ul + 1568 (p 21 +p 24 ) (B35) 

where ^ is the control input to the inboard elevon and (P2i + P24^ is t * ie SAS 
input to the inboard elevon. 


Inboard actuator 



Figure B2. Inboard/ Outboard Elevon Actuator 
Block Diagram 
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The governing differential equations for the outboard elevon system fall 
into two separate cases. The first case deals with the standard system with 
only an additional control input, Ug, to the outboard actuator. The governing 

differential equation for this case is readily determined from Figure B2 as: 

6 2 = -67.8 fc 2 - 1149 6 2 + 1149 6 + 1149 u 2 (B36) 

The transfer function for the second case is shown in Figure B3. This 
system is designed to allow independent operation of the inboard and outboard 
elevon system by canceling the mechanical linkage input of 6^ to the outboard 

actuator. The differential equation for 6 2 (t) is found to be: 

<5 2 = -67. 8 8 2 - 1149 § 2 + 1149 p 23 + 1149 u 2 

p 23 = P 22 (B37) 

P 22 = ~ 50 - 5 P 22 " 1568 p 23 + 1568 ^ P 21 + P 24^ 



Figure B3. Electronic Cancellation of Mechanical Linkage 
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Canard and mass -reaction device state equations. - The canard actuator 
state equation for the optimal control studies was defined to be: 


6 3 = 30 S 3 + 30 u 3 

where u 3 is the control input into the canard actuator. 

The dynamics for the mass-reaction device were defined as: 

6 =-2 6 - 4 6 - 4 u 

m m m m 

where u m is the control input to the actuator. 


(B38) 


(B39) 


Baseline pitch SAS state equations . - The transfer function for the base- 
line pitch SAS is shown in Figure B4. For the frequency range of interest in 
this study, one can neglect the rate gyro and servo dynamics. For F/Cs 1, 
1A, and 3, one has the following differential equation for the SAS output, p • 

16 

p 2 i - °- 49 5 e « 0 + 3 - 94 V s + °- 49 6 e I 0 j ' ( W ' r j + 2 ) 

j=l 

-3. 94 p 21 (B40) 

where 


0 , 


r . 
1 


scheduled gain 
til 

slope of j mode shape at the rate gyro location 
til 

j generalized coordinate in eq. (B25) 


For the above flight conditions, p 24 = 0. 

For F/C 2, the SAS output has an additional term, p 24> corresponding to 
the lagged pitch rate path in Figure B4. The governing state equations for 
P 2 4 are: 16 

P25 = 4 - 006 ^o +4 ' 006 b o + 4 - 006 I V fr 3+2 + V j+2 ) 

3=1 

- 10. 49 p 25 - 0. 8012 p 24 (B41) 


P 24 


= P 


25 


where 0 , 0/and r. are defined as before, 
o J 1 
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Figure B4. YF-12A Baseline SAS Block Diagram 
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Total state vector and state equations . - Using the new state variables and 
associated state equations defined in the preceding sections, eq. (B25) can be 
readily transformed into the required form using standard techniques. Thus, 
the form of the state equations becomes: 




x = Fx + G^u + G 2 r) 

(B42) 

where x^ = 

Lw, 0, r 2 , . . 

** r 18’ r 2’ ‘ ‘ * ' r 18 > P 20 J W g' 



5 r * 5 2’ 

5 2' p 22' p 23’ 6 3' ^m' 6 m^ 


T 

u = 

Lu 1 , u 2 , u 3 . 

u J 
m 


r] = 

white noise 




One notes that in the above state vector, states p 21 . 


p 24 , and p 25 are 


missing. These states are associated with the SAS and will be added after the 
feedback states have been transformed into accelerations and rates. Further 
observation shows that both the canard and mass -reaction state equations 
appear in eq. (B42). For any of the problems considered in this study, at 
most one of these devices will be used, and hence eq. (B42) must be modified 
to reflect this at the time of its use. 


Transformation of the feedback states. - It is desired that the feedback 
states be velocities and accelerations rather than displacements and veloci- 
ties. The following transformation will accomplish this objective but is 
different for the case of a mass-reaction control instead of a canard control. 
Thus, it is necessary to state two transformations, and hence we shall define 
two new state vectors and the associated state equations. 


Let 


[x (1) ] = ie Q , w q , r g , . . . , r lg , r 3 , . . . , r lg , p^ . 


20 ’ 


w r 5 1’ 5 r 5 2' 5 2' p 22’ p 23' 5 3 J 


(B43) 

be the state vector without the mass-reaction control satisfying the equation 

(B44 ) 


x (1) = F (1) x (1) + G. (1) u (1) + q 


where F^, G^\ are same as F, G^, G 2 with columns and rows 

1 and 2 interchanged along with the removal of the states and equations asso- 

(1) T 

ciated with the mass-reaction device. The [u v ] term in this case is just 
Lu 1 , u 2 , u 3 J. 
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Similarly, let 

(2) T 

[x ] = L0 o , w o , r 3 , .... r 18 , r 3 , r 


18' P l’ • * ‘ ’ p 20' 


w g , 6 1> 6 1 , 5 2 , 6 2 , p 23 , p 24 , 6 m , 6 m J 


(B45) 


be the state vector for the system without the canard but including the mass- 
reaction control satisfying the equation 


( 2 ) = f ( 2 ) x ( 2 ) + g i ( 2 ) + g 2 ( 2), 1 


(B46) 


( 2 ) ( 2 ) ( 2 ) 

where F' , , Gg are the same as F, G^ G 2 withrows and columns 

1 and 2 interchanged along with the removal of the canard states and equations. 

( 2 ) ^ 

The [u '] term thus becomes 

[u (2) ] = Lu r u 2 , u m J 

The transformation for the system with a canard controller is developed 
by first partitioning the state vector x^ as follows: 

x (1) = (X, Y, Z) (B47) 

where 

X = 6 

o 

Y = Lw q , r 3 , . . . , r lg , r 3 , . . . , r lg J 

Z = LPj, p 20 , Wg, 6^ 6 2 , 6 2 , p 22> p 2g , 6 g J (B48) 

Thus, the governing differential equation for this system can be written as: 


x' 


1 

* 

F 12 

< 

CO 

r 

> ■ 

F 2 1 

F 22 

F 23 

z 

\ 


F 31 

F 32 

oo 

CO 


V 


"0 


f 

0 

) 

+ 

0 


° 

z 

\ 1 

1 

G 1 3 


° 23 , 


(B49) 
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Now 


I O O 


*21 *22 F 2 V Y 


O O I 


TF (1) < Y 


(B50) 


TF (1) < Y 


= TF^ 1 ' F^ TF^ 1 * 


y U F <y g 13 {u ( 1) }+ wy g 23 U 


i<« = F< X > X (1) + G™ u (1) + G ? (1) n 
new new new 1_„„ 

new new 


r(l) . 

^new ^ = ^o’ w o' r 3’ '•*’ r 18' r 3' ' * * ' r 18’ p l* * * * * p 20’ 

W g' ®1’ ' 6 2 j ®2 j P 22’ p 23’ ^3 -1 

and F G, , G„^^ are defined by eqs. (B51) and (B52). 
new 1 2 J ^ 

new new 

The transformation for the system with a mass-reaction controller is 

(2) 

achieved by first partitioning the state vector x as follows, where: 
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X = 6 


- Lw 


o' 3' 


18' 


r 18 J 



p 20’ W g* 6 1’ 6 1' ®2* 5 2’ P 22’ p 23' ®m* 


(B54) 


One notes that when the mass -reaction control is used in the system, the 
matrix G j 2 ^ contains elements in the first 18 rows that do not appear in . 

Thus, we now write: 


U l' U 2 

u 

m 


1 


— 

— 


- | 

— 

O 

r (2) 

G llm 

= 

A(2 ) 1 

G 1 i 

g< 2 > 

lm 

O 

G 12m 


1 


'(2) 

_13 

r (2) 

G 13m 


1 

1 





1 

— 


(B55) 


Thus, the state equation becomes: 


' X ' 


r F (2) 

F n 

f ( 2 ) 

F 12 

*1 

i 


r X' 


|' G llm 


p) 

\ 


' O ' 


Y 

>- 

p(2) 
F 2 1 

f ( 2 ) 

F 22 

f ( 2 ) 

F 23 


Y 


A(2) 1 

G 1 1 G 12m 

< 

< 

U 2 

► 

> + < 

o 

) T) (B56) 



JL2) 

.31 

f ( 2 ) 

F 32 

f ( 2 ) 

F 33_ 




1 G 13m_ 

' 

u 

m J 


r (2) 
v u 23 J 



One now solves for 



and obtains 


x' 


"I 

O 

° 1 

-1 

r X ' 


‘I 

O 

O ' 

-1 

' O 

Y 

> = 

f (2) 

F 21 

f (2) 

F 22 

^ CM 

( 

Y 

> + 

f (2) 
F 2 1 

f (2) 

*22 

f (2) 

F 23 


G 12m 

k Zj 


O 

O 

. J 


, Z , 


O 

O 

I - 


-O _ 


(B57) 
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Let 


TF 


( 2 ) _ 


I 

f (2) 
F 2 1 

O 


O 


O 


,( 2 ) 

22 


,< 2 ) 

23 


O 


(B58) 


Thus: 



-1 


(x) 


Y 


v. 


z 


J 


{ G 1 2 m ) U m 


1 ° ) 


r~ — 


( \ 



G llm| 


' 

G 12m 

_ 


G 13m ; 


'12m, 




m 


1° J 


u 


m 


(B59) 


The above state equation can be simplified by adding u m as a state and 
letting u^ be a control. Thus, one obtains the following state equation: 
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(B60) 


where 


F = [TF (2) ] [f (2) ] [tF (2) ] 

Thus, the new state equation becomes: 


^(2) = f (2) x ( 2) 
new new new 

, r (2) (2 ) G _ 
1 new u new 2 new 

(B61) 

where 



, (2) .T , a 

x new " ^ 6 o > W o J r 3' * * * ' 

• • • • 

r 18' r 3' ’ * • J r 18' P 1 J **** P 20 J 


W g» ®2' 

5 2' P 22' P 2 3' 6 m’ U nJ 

(B62) 

and 



“new* = L“i- “2- “m 


(B63) 


Once the system has been transformed, the SAS states which are func- 
tions of 0 q , m and m can easily be appended. The remaining task is to derive 

the response equations. 
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Description of the response equations. - The response equation used in 
the optimal control analysis is of the form 

r = Hx + Du (B64) 

where x is the state vector and u is the controller. 


The desired responses, used in the study, are the vertical accelerations 
in g at FS 248, FS 738, and FS 1236 at WS32, along with the control vector. 
The vertical accelerations at the above points are given by: 

1 J.*_ - « . - *« . Y_ 0. r. + 2 |j = 4, 8, 16 (B65) 

i=l j > 


n 


386 


\ w - 

( ° 


v e + x. e + 

O 3 


where 

— i- 

x. = distance in inches from the j station to the reference point, 
^ FS 900, and is negative for points forward of FS 900 

i.U 

0. = deflection at station j in the i elastic mode 


The desired response vector is of the form: 

A -p 

(r) = Ln , n , n , u , u 9 , u~ or u J (B66) 

4 Z 8 z 16 l4o m 

Hence, from the definition of r and n , one can form the matrix equation 

3 

A. 

r = Hx+Du (B67) 


Optimal practicalization technique. - In previous LAMS applications the 
feedbacks specified in the optimal controllers included not just the mode feed- 
backs (including the rigid-body mode) but also the actuator states and the wind 
model states. Since the wind model states are not truly measurable, it was 
desirable to eliminate these feedbacks from the optimal controllers. Further- 
more, it was desirable to generate optimal controllers which also did not con- 
tain some of the measurable feedbacks such as the rigid-body modes. To 
provide this increased flexibility in the optimization program, an additional 
optimization procedure was applied. This procedure is referred to as the 
practicalization procedure. It is used to drive the gains associated with the 
undesirable feedbacks to zero. A suboptimal controller is achieved by inte- 
grating a parameter, X, from 1 to 0 in the equation: 

K(M = K 1 (X) + \K 2 (B68) 

while minimizing the cost J*. In the term [K(M], K^X) includes the gains on 
those feedbacks to be retained, and K 2 includes those gains on the feedbacks 
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to be discarded. This is done by computing dK^/dX using the implicit function 

theorem and numerical integration techniques. This practicalization scheme 
was used successfully in the Phase I study and is considered as a valuable 
tool. 


Integration of synthesis techniques. -- Application of quadratic optimiza- 
tion techniques yields a controller which provides the required allocation of 
control effort among the force producers. To reiterate, however, the dis- 
advantage of the technique is that the feedbacks are specified in terms of the 
state variables, rather than in terms of practical sensor outputs (e. g. , 
accelerometer or rate gyro). The remaining task, then, is to make the 
transformation from state variable feedbacks to feedbacks from sensor out- 
puts. The sensor output signals are simply linear combinations of the state 
variables and each combination is determined by the sensor location on the 
vehicle. 

There are two schemes which can be used to make the transformation 
from state variable feedbacks to feedbacks from practical sensors. The first 
scheme consists of the following four steps: 

1. Eliminate force producer crossfeed signals. For example, 
eliminate outboard elevon deflection and rate signals into the 
inboard elevon controller and vice versa. This is achieved by 
a simple matrix transformation which results in establishing 
filtering on each force producer input. 

2. Eliminate insignificant state variable feedbacks. This is 
achieved by computing frequency responses of each feedback 
variable for a force producer input. The frequency responses 
are then compared with each other and with unity to determine 
their relative significance. This analysis is illustrated in 
Figure B5. The insignificant feedbacks are eliminated and 
stability is checked with the remaining feedbacks. 

3. Reorganize remaining state variable feedbacks into groups so 
that each group is a linear combination of variables resembling 
a rate gyro or accelerometer output. 

4. For each group in Step 3, find a sensor combination which will 
yield similar coefficients on the state variables. This is 
largely a trial and error procedure. One can always use as 
many sensors as there are state variables in each group and 
compute an ideal blend of the sensors. In most cases, however, 
this is not practical, especially when tolerance to parameter 
variations is taken into account. 

The second scheme for transforming from state variable feedbacks to 
practical sensor feedbacks is to use the frequency response synthesis tech- 
nique. In this scheme, the attenuation factors, T\, required in the frequency 
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6^ = force producer deflection 
ri . = state variable 



1+[K 1 ( 


^♦k 2 <? 
1 c °1 
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Figure B5. Elimination of State Variable Feedbacks 


response technique are computed using the optimal control system. This, in 
essence, establishes the desired allocation of control effort which is lacking 
in the frequency response technique. Having computed the TVs, the frequency 

response technique is applied to synthesize a control system. The force 
producer crossfeed signals in the optimal controller must be eliminated prior 
to computing the IVs. The method used to compute the TV's using the 

optimal controller is outlined in the following paragraphs. 

In order to use the frequency response technique in the synthesis of 
optimal controller, it is necessary to make the appropriate transformations 



a form compatible with the frequency response technique as described in this 
appendix. Essentially, this means finding the appropriate T's at the fre- 
quencies of interest, such as the resonant peaks. 

The inboard/ outboard elevon controller of Figure B6 shows the corres- 
ponding terms of the optimal controller and the frequency response techniques 
where : 
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• Hj and Hg are the transfer functions obtained from the 
optimization. 

• Gj and Gg are the filters of the frequency response technique 
with their respective sensors, z t and z 9 . 

SI S ct 

The equations presented below are those for the 6^ controller, but simi- 
lar equations can be derived for the 6 ^ loop. 

From Figure B6, it can be seen that 


H 


1 



6 2 CL 


or by solving for G^(s) 


K^s) 


1 - K i G i (8 >i* 


si 


6 2 CL 


(B69) 


KjG^s) = 


Z sl /6 1 


6 2 CL 


r 

i 1 

i 

h - 4 

+ 

1“ 

Z Sl(> 


*sV 8 l] 

6 2 cl _ 


(B70) 


There is also a stability constraint on the above equation as evident from 
Figure B6. The stability of the 6^ loop is determined by the following equa- 
tions: 


KjG^s) 



6 2 CL 



(B71) 


The r of the 6^ controller can be determined at the frequencies of interest 
from eqs. (B4) and (B70): 



(B72) 


From Eq. (B72) it can be seen that it is possible to compute r^, knowing 
and the frequency responses of acceleration to gust and surface inputs at 
the vehicle station of interest. Knowing r^, it is then possible by using the 
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frequency response technique to determine the necessary filter and sensor 
combination for the 6^ loop. Substituting G.^ and Z g ^ into the aircraft 

equations, one may then synthesize the 5 ^ control loop by a similar proce- 
dure. Therefore, this method allows one to synthesize an optimal controller 
on a loop-by-loop basis. 


Synthesis of Optimal Controllers 

An attempt was made in the study to use both of the above schemes for 
transforming the state variable feedbacks to practical sensor outputs. Unfor- 
tunately, the optimal controller which was generated was not a reasonable 
system in terms of rigid-body handling qualities and stability constraints. 
There was not sufficient time nor justification in the study to warrant rede- 
fining the optimal controller. What is required is a handling quality constraint 
in the optimization cost function. The optimal controllers which were defined 
are described in the following paragraphs. None of these made use of a 
handling quality constraint, however. Instead, attempts were made to obtain 
reasonable handling qualities by trial and error adjustment of the quadratic 
weights in the cost function. This was not successful. 

Optimal quadratic controllers were computed for various combinations of 
inboard and outboard elevons, and canard or mass-reaction force producers. 
These controllers were computed using selected quadratic weights at the high 
subsonic, heavyweight flight condition. Optimal controllers were also com- 
puted for a single set of quadratic weights for all four flight conditions using 
the inboard and outboard elevons. 


Table B1 shows the ideal gains obtained with the inboard and outboard 
elevons at each of the four conditions. As is evident from this table, only 
five structural modes were included in F/Cs 1, 1A, and 3, but six modes 
were included at F/C 2. The p lg , p 2Q , and w g terms represent wind model 

states which are usually eliminated by applying optimal practicalization pro- 
cedure. The 6 2 are inboard elevon position and outboard elevon posi- 


tion, respectively. Hence, each controller contains crossfeeds between the 
inboard and outboard elevons. These crossfeeds can be eliminated by a 
simple transformation. The pg^, $ 24 ’ an< * ^21 * erms represent states of the 


baseline SAS which are measurable. They also can be eliminated if so desired 
by the optimal practicalization procedure. 


Table B2 shows the rms accelerations of these controllers at the four 
flight conditions. These results show substantial ride improvement is pos- 
sible with these controllers at all four conditions. Figures B7, B8, and B9 
show the PSD plots for the high subsonic, heavyweight flight condition. These 
responses show considerable rigid-body load alleviation in addition to good 
mode suppression performance. 
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TABLE B2. - RMS ACCELERATION PERFORMANCE WITH AN 
OPTIMAL CONTROL SYSTEM USING INBOARD AND 
OUTBOARD ELEVONS (QUADRATIC WEIGHTS = 

2 , 1 , 1 , 1 , 1 ) 


Fuselage 

station 

F/C 1, 

high subsonic, 
heavyweight 
condition 

F/C 1A, 
high subsonic, 
lightweight 
condition 

F/C 2, 
cruiBe, 
heavyweight 
condition 

F/C 
low su 
lightw 
cond 

— 

: 3, 
bsonic 
eight 
tion 

SAS 

only 

LAMS 

SAS 

only 

LAMS 

SAS 

only 

LAMS 

SAS 

only 

LAMS 

248 (pilot) 

0. 0204g 

0.008 g 

0. 022 7g 

0. 0083 g 

0. 0095g 

0. 0036g 

0. 0085g 

0. 0033 g 

738 

0. 0131 

0. 0052 

0. 0181 

0. 0065 

0. 0054 

0. 0035 

0. 0092 

0. 0035 

1236 (aft fuselage) 

0. 0211 

0. 0060 

0. 0031 

0. 0080 

0. 0108 

0. 0038 

0. 0121 

0. 0052 


At F/C 1 (i. e. , the high subsonic, heavyweight condition) the optimal 
practicalization procedure was used to eliminate various feedbacks from the 
optimal control system. Of particular interest is the case where only the 
wind model states are eliminated. 

Table B3 compares optimal gains with and without the wind states (p^g, 

P 2 Q» w ). Figures B7, B8 and B9 show the change in the mean square 

acceleration performance resulting from elimination of the wind states. Also 
shown are gains obtained with elimination of the baseline SAS signal (p 91 ) and 

the rigid-body feedbacks (w , 0 Q ). To observe the effect of removing the wind 

states (which are not measurable by a sensor) it is best to examine the change 
in the closed-loop poles. Table B4 shows the roots of the rigid-body and 
structural modes for the optimal system with and without the wind states. 

The rms accelerations increased only slightly when the wind states were 
removed. But. as Table B4 shows, the short-period and the first structural 
mode roots changed significantly. The values of the roots shown indicate that 
the handling qualities of the airplane would be significantly affected by this 
controller. These results tend to substantiate the hypothesis that the cost 
function should be modified to include a handling qualities constraint such as 
model following. Additional practicalization of this optimal controller by 
eliminating the baseline SAS feedbacks to the system, and the rigid-body 
feedbacks failed to improve the character of these roots. 

An effort was made to replace the idealized feedbacks with feedbacks 
from practical sensors for the optimal system without wind states. Begin- 
ning with the feedbacks shown in Table B1 for the system without wind states, 
the first step is to eliminate the actuator crossfeeds between the inboard 
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Figure B7. Mean Square Acceleration at FS 248 for a 1-ft/sec Turbulence 
Input With and Without Optimal Controllers -- F/C 1 
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Figure B8. Mean Square Acceleration at FS 738 for a 1 -ft/ sec Turbulence 
Input With and Without Optimal Controllers -- F/C 1 
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Figure B9. Mean Square Acceleration at FS 1236 for a 1-ft/sec Turbulence 
Input With and Without Optimal Controllers -- F/C 1 
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TABLE B4. - COMPARISON OF ROOTS FOR AN OPTIMAL 
SYSTEM WITH AND WITHOUT WIND STATES 


Root 

description 

Open-loop roots 

Optimal system 
with wind states, 
closed-loop roots 

Optimal system 
without wind states, 
closed- loop roots 

Damping 

ratio 

Frequency 
(rad/ sec) 

Damping 

ratio 

Frequency 
(rad/ sec) 

Damping 

ratio 

Frequency 
(rad / sec) 

Higid body 

0. 62 

3. 77 

0. 607 

3. 30 

0. 910 

2.05 

Mode 1 

0. 08 

12.43 

0.224 

9. 12 

0. 482 

5. 97 

Mode 2 

0. 054 

16. 1 

0. 851 

10. 88 

0. 648 

12. 87 

Mode 3 

0. 061 

28. 7 

0. 527 

36. 85 

0. 624 

24.21 

Mode 4 

0. 040 

44. 47 

0. 941 

40.28 

0.312 

44. 60 

Mode 5 

0. 060 

50. 07 

0. 044 

47. 82 

0. 079 

48. 57 


elevon and outboard elevons. This step, after some additional simplification, 
results in the system shown in Figure BIO. The input signals z g ^ and z g 2 

are defined as: 


z 1 = -0. 00291 z. + 1 0- (B73) 

si 1 1 

where 

z, = z - 0. 302 r Q + 0. 143 r. - 0. 088 r, + 0. 0335 r~ - 0. 134 v n 
1 o o 4 o b i 

9 X = 0. 515 9 q + 0. 0 146 r 3 - 0. 0066 r 4 - 0. 0215 r g - 0. Oil r g + 0. 006 r ? 

and 

z s2 = 0. 00306 z g + 9 2 (B74) 

where 

z 2 = z q - 0. 046 r 3 - 0. 145 r 4 + 0. 220 r g + 0. 063 r g - 0. 029 r ? 

0 2 = -0. 743 0 Q - 0. 0028 r g - 0. 0085 r 4 - 0. 0168 r g + 0. 006 r g - 0. 001 r ? 
The optimal feedbacks from the baseline SAS output, p_ , appear to 

u 1 

negate the effects of the baseline SAS, at least on the inboard elevon. Sub- 
sequent analysis indicated, however, that these two p 21 feedbacks in the 
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optimal system could be eliminated without a significant change in the ride 
quality performance. It is of interest to note that this optimal system con- 
tains compensation very similar to the compensation derived by frequency 
response analysis for the inboard/ outboard elevon control system. The 
0. 2 -sec first-order lag is predictable also from the application of North 
American's ILAF concept (see ref. 3). 

As discussed previously, further attempts to reduce this controller to a 
system using practical sensors was not successful. The nature of the feasi- 
bility study did not warrant further analysis of this problem in the study. 

Several optimal controllers were defined using either a canard or mass- 
reaction configuration plus the inboard and outboard elevons. However, in all 
cases the same problem was encountered. Reasonable handling qualities 
could not be assured. Regardless of the weighting factors selected, the rigid- 
body and first-mode frequencies were reduced to unacceptable values. Root- 
mean-square performance with either the mass-reaction configuration or the 
canard configuration was substantially better than that obtained with just the 
inboard and outboard elevons. 
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